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TRANSITION FROM GAUSSIAN TO NON-GAUSSIAN 
FLUCTUATIONS FOR MEAN-FIELD DIFFUSIONS IN SPATIAL 

INTERACTION 

ERIC LUQON^ AND WILHELM STANNAT^’® 


Abstract. We consider a system of N disordered mean-field interacting diffusions 
within spatial constraints : each particle 6i is attached to one site Xi of a periodic lattice 
and the interaction between particles 6i and 9j decreases as \xi - Xj\~°‘ for a e [0,1). In a 
previous work [^, it was shown that the empirical measure of the particles converges in 
large population to the solution of a nonlinear partial differential equation of McKean- 
Vlasov type. The purpose of the present paper is to study the fluctuations associated to 
this convergence. We exhibit in particular a phase transition in the scaling and in the 
nature of the fluctuations: when a e [0, ^), the fluctuations are Gaussian, governed by 
a linear SPDE, with scaling '/N whereas the fluctuations are deterministic with scaling 
N^^°‘ in the case a e (L, 1). 

2010 Mathematics Subject Classification: 60F05, 60G57, 60H15, 82C20, 92B25. 


1. Introduction 


The aim of the paper is to study the large population fluctuations of disordered mean- 
field interacting diffusions within spatial interaction. A general instance of the model may 
be given by the following system of N coupled stochastic differential equations in X := 

(m ^ 1): 


dOi^t = c{6i^t,^i)dt + 


rn T. r(«.,i Xj) dt + dBi^t, 

i*i 


i e Aat, t ^ [0, T] 
( 1 . 1 ) 


where T > 0 is a fixed but arbitrary time horizon. In p.ll) . c{6i^t:^i) models the local 
dynamics of the particle 6i, T{6i^t,i^i,9j,t,^j) governs the interaction between particles 6i 
and 6j and is a collection of independent standard Brownian motions modeling 

thermal noise in the system. Both local dynamics and interaction are perturbed by an 
independent random environment that is a sequence of i.i.d. random variables (a;*)* i in 
Y := R” (n ^ 1) modeling some local inhomogeneity for each particle. 

The novelty of dni) is that some geometry is imposed on the interactions: the particles 
9i in dLH) are regularly positioned on a periodic lattice (of dimension 1 for simplicity) 
and the interaction between two particles depend on the distance between them. Namely, 
define Aw := {-N,... ,N} with -N and N being identified, with cardinal |Aw| = 2A and 
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suppose that for all i e A^r, the particle 9i is at the fixed position Xj := ^ e §n, where 
E>n ■= {xi ■= i e Atv} is a subset of the one-dimensional circle S := M/Z. Making the 
obvious identification between x e [-^,1] and its equivalence class x e S, the euclidean 
norm I’l on induces a distance on § by: 

Vx = x[l],y = ?/[1] e§, d{x,y) = m.m{\x-y\,l-\x-y\). (1.2) 

The interaction in (HI) is supposed to decrease polynomially with the distance between 
particles through the spatial weight 

T(x,y) :=(i(x,y)"", x,yeS, (1.3) 

where a is a parameter in [0,1). 

Remark 1.1. We will often drop the notation j i in sums as in dni), sinee, with no 
loss of generality, one can define T(x,x) ;= 0 for all x. 

1.1. Weakly interacting diffusions. 

Mean-field models in physics and neuroscience. In the particular case where a = 0 and 
Ui = 0, the geometry and the disorder in (jl.ip become irrelevant and we retrieve the usual 
class of weakly interacting diffusions (studied since McKean, Oelschlager and Sznitman [311 
UBIIMIEH]). Such systems are encountered in the context of statistical physics and biology 
(synchronization of large populations of individuals, collective behavior of social insects, 
emergence of synchrony in neural networks PSollIl]) and for particle approximations of 
nonlinear partial differential equations (see e.g. [inii29] and references therein). Usual 
questions for this class of models concern their large population behavior (propagation of 
chaos, existence of a continuous limit for the empirical measure of the particles [MiiiaET], 
fluctuations around this limit [MlEHlET] and large deviations m) as well as their long¬ 
time behavior P|29]. 

A recent interest for similar mean-fields diffusions has been developed for the modeling of 
the spiking activity of neurons in a noisy environment (e.g. Hodgkin-Huxley and FitzHugh- 
Nagumo oscillators [6l[T5l[l0]). In this context, 6i represents the electrical activity of one 
single neuron, F captures synaptic connections between neurons and the disorder Wj models 
an inhomogeneous behavior between inhibition and excitation. We refer to [6] for more 
details on application of these models to neuroscience. 

Another illustration of weakly interacting diffusions concerns statistical physics and 
models of synchronization of oscillators. In particular, the Kuramoto model (see e.g. 
[H1271 dal ET]) describes the synchronization of rotators 9i in M/27rZ with inhomogeneous 
frequencies Wj: 

K ^ 

‘i9i{t) = cjjdt H-^ sin {9j - 9i) dt + adBfit), tfid, i = 1,... ,N. (1-4) 

^ j=i 

The system (11.41) is well-known to exhibit a phase transition between incoherence to syn¬ 
chrony as the interaction strength K > 0 increases. We refer to the mentioned references 
for further details on the dynamical properties of (|1.4D . 

Diffusions in spatial interaction. The motivation of going beyond pure mean-field inter¬ 
actions comes from the biological observation that neurons do not interact in a mean-field 
way (see mm and references therein). There has been recently a growing interest in 
models closer to the topology of real neuronal networks [21( [35l 141] . Even though the 
analysis of such models seems to be difficult in general, it is quite natural to expect that 
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properties valid in the pure mean-field case (the first of them being the existence of a 
continuous limit in an infinite population) still hold for perturbations of the mean-field 
case, namely for systems where the interactions are not strictly identical, but where the 
number of connections is sufficiently large to ensure some self-averaging as the population 
size increases. 

In this perspective, the main motivation of (ll.ljl comes from works of Rogers and Wille 
[36] and of Gupta, Potters and Ruffo m (see also mn where a spatial version of the 
Kuramoto model is introduced; 

K ^ 1 

d9i^t = dt + — ^ sin {6j t - Oi^t) ■, - dt + a dRj t ^ 0, z = 1,..., iV. (1.5) 

^ pi 

In (|1.5I1 (and more generally in (ll.ip i. the particles are still interacting on the complete 
graph but the strength of interaction decreases polynomially with the distance between 
particles. 

System (|1.5p has to be related to analogous models of statistical physics with long-range 
interactions (e.g. the Ising model or XY-model with interaction in r““, see [3] and refer¬ 
ences therein). The influence of the spatial decay in (11.51) on the synchronization properties 
of the system (in particular the existence of critical exponents and finite size effects) does 
not seem to be completely understood so far (see the aforementioned references for further 
details). 

Empirical measure and McKean-Vlasov limit. Note that in (jl.ip . e X ;= R”^, but one 
should also think to the case of X being a compact domain of R”^, see for example X = 
R/27rZ in the Kuramoto case (II.5p . For any vectors u and u in X (or Y), |ri| stands for the 
euclidean norm of u and u ■ u is the scalar product between u and v. 

We endow (jl.ip with an initial condition: the particles (0i,o)i=i,...Af are independent 
and identically distributed according to some law ( on X. The disorder (wi)i=i,...,Ar is a 
sequence of independent random variables, identically distributed according to some law 
H on Y, independently from (6'i,o)ieAjv and (Rj)ieAjv 

All the statistical information of (HH) is contained in the empirical measure of the 
particles, disorder and positions (that is a random process in the set of probabilities on 
Xx Yx§): 

E te[0,r]. (1.6) 

The object of a previous work |28j was to show that, under mild assumptions on the model 
(see [28], Theorem 2.18), the empirical measure converges, as W ->■ oo, to the unique 
solution u of the following weak McKean-Vlasov equation: 

//<*"' = //d-'o+X f + J r(-,0,£i))T(',x)zzs(d0, dia, dx)|jdi^sds, 

(1.7) 

for any regular test function {6,uj,x) f {9, lj,x) and where the initial condition is 

uo( d9, dw, dx) ;= C{d9)fi{ du) dx. (1.8) 

The main difficulty in the analysis of dEZl) is that it involves the singular kernel T so 
that even the well-posedness of such equation is unclear. We refer to [28] for more details. 
Equation (jl.7|) is a spatial generalization of the standard McKean-Vlasov equation in the 
pure mean-field case (see e.g. [331 dl ETj). 
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Fluctuations around the mean-field limit. The purpose of this work is to address the ques¬ 
tion of the fluctuations of the empirical measure visf (II.6p around its limit n (jl.7p . We are 
in particular interested in the influence of the geometry on the speed of convergence of visf 
towards n. 


Definition 1.2. Define the fluctuation process 

'nN,t ■= aN {TyN,t - , N^l, te[0,r] (1.9) 

where on is an appropriate renormalizing factor. 


The first contribution of the paper is to exhibit a phase transition in the size of the 
fluctuation of m- the correct scaling aw in (|1.9p depends on the value of the spatial 
parameter a e [0,1) in (11.31) with respect to the critical value 


Definition 1.3 (Fluctuation renormalization). Fix 0 ^ a < 1 and define for N 1 


( 1 . 10 ) 


Let us show intuitively that Definition 11.31 provides the correct scaling for (jl.9p . The 
convergence of towards n is due to the competition of two effects: the convergence 
of the empirical distribution of the initial condition Oi^ and the Brownian motions Bi 
(which scales typically as \/]V) and the convergence with respect to the spatial variable 
Xi. To fix ideas, set c e 0 and T e 1 in (jl.ip . In this case, everything boils down to the 
approximation of the integral /g dx by the Riemann sum T,f=i • A simple 

estimate (see Lemma 13.41 below) shows that the rate of this last convergence is exactly 
Definition 11.31 simply chooses the predominant scaling in both cases. 

This intuition also suggests that when a < |, the randomness prevails and one should 
obtain Gaussian fluctuations as oo, whereas when a > | the randomness disappears 
under the scaling and one should obtain a deterministic limit for r/w- The main 

result of the paper is precisely to make this intuition rigorous: we show that the fluctua¬ 
tion process r/w converges to the unique solution of a linear stochastic partial differential 
equation when a < ^ (see Theorem EH) and that rj]^ has a deterministic limit in the 
supercritical case a > ^ (see Theorem 12.Sp . 

What makes the analysis difficult here is the singularity of the spatial kernel T in (II.ip . 
An important aspect of the paper is the introduction of an auxiliary weighted fluctuation 
process V-n (in addition to r/w) that is necessary to capture the spatial variations of the 
system and to cope with the singularity of the weight T. We refer to Section E] for more 
details. 


On the quenched fluetuations. The main results (Theorems E21 and EH are averaged with 
respect to the disorder (wj)*. Although we did not go in this direction for the simplicity 
of exposition, analogous results also hold in the quenched set-up, that is when we only 
integrate in (11.111 w.r.t. the Brownian noise and the initial condition, and not w.r.t. the 
disorder. This has been carried out in a previous paper m in the non spatial case. We 
let the reader adapt a similar strategy to the present situation. 

Existing literature. The use of weighted empirical processes such as JIn (see Section EH 
in the context of interacting particle systems is reminiscent of previous works. One should 
mention in particular the articles of Kurtz and Xiong [25( 126] on particle approximations 
for nonlinear SPDEs. 
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The present paper uses Hilbertian techniques already introduced by Fernandez and 
Meleard |16] (see also |38[ 127] ) who proved a similar central limit theorem for weakly 
interacting diffusions without spatial geometry (that is a = 0 in the framework of (|l.ip i. 
It is shown in [16] that the fluctuations are governed by a linear SPDE and that the 
convergence holds in some appropriate weighted Sobolev space of distributions. One can 
see the first result of the paper (Theorem [22]) as a generalization of the result of Fernandez 
and Meleard to the spatial case: when a < ^, the spatial damping on the interactions in 
dEU) is not strong enough to have an effect on the behavior of the fluctuations of the 
system and the fluctuations remain Gaussian. 

The supercritical case (Theorem (221) may also be related to a class of models previously 
studied in the literature, that is moderately interacting diffusions (see [23l|32]). This class 
of models also exhibits deterministic fluctuations (see [23] p. 755), but one should point 
out that the precise scaling oat is not explicitly known in this case (see [23], Remark 3.15). 


Comments and perspective. The main conclusion of Theorem 12.81 is that, in the case a > |, 
the leading term in the asymptotic expansion of the empirical measure around u is 
deterministic, of order . A natural question would be to ask about the existence 
and the nature of the next term in this expansion. Concerning the dependence in the 
spatial variable, one easily sees that the term following in the asymptotic expansion 
of is of order ■^, which is in any case smaller than the Gaussian scaling 

Consequently, one should expect the next term in the expansion of UAr to be Gaussian, of 
order The precise form of this term remains unclear, though. Note that a similar 

analysis has been made by Oelschlager in [34] concerning moderately interacting diffusions. 

Another natural question would be to ask what happens at the critical case a = | (that 
is when the spatial and Gaussian fluctuations are exactly of the same order). Although it 
is natural to think that the correct scaling is exactly \/N, the present work only provides 
partial answers. A closer look to the proof below shows that the scaling is at least 
(see Section 12.4p . To derive the correct scaling and limit for the fluctuations in this case 
seems to require alternative techniques. 

The behavior as N ->■ oo of systems similar to (II.ip in the case a > 1 is also of interest 
and is the object of an ongoing work. The existence of a continuous limit at the level of 
the law of large number is doubtful in this case. 

It is likely that the results presented here should be generalized to other models of 
interacting diffusions, especially to systems with random inhomogeneous connectivities 
(in the spirit of [HIT]) which are of particular interest in the context of neuroscience [6]. 
This will be the object of a future work. 


Outline of the paper. In Section |21 we specify the assumptions on the model and state 
the main results (Theorem 12.71 and 12.8p . Section |3| is devoted to prove tightness of the 
fluctuation process in an appropriate space of distributions. The identification of the limits 
is done in Section Hi Section [5] contains the proofs of technical propositions. 


2. Main results 

In the rest of the paper, r = (0,uj,x) stands for an element of X x Y x S. In particular, 
for any s e [0,T], i e Atv, Ti^s = We use also the duality notation 

(A,/):= f f{6,uj,x)X{d6,du,dx), 

JXxYxS 


( 2 . 1 ) 
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where A is a measure on X x Y x Z and {9, uj, x) f{6, co, x) is a test function on X x Y x S. 
We also use functions of two variables defined on (X x Y x S)^ and introduce notations for 
the corresponding partial integrals; for any such function (r, f) g{T,f), define 


[ 9 , A ^(r) = [ 5 , A ^{9,u},x) := 
[ 9 , = [ 9 , X]^{0,oj,x) := 


g{9,u},x,9,Cij,x)X{d9, duj, dx), (2-2) 

g{9,u},x,9,Cij,x)X{d9, duj, dx). (2-3) 

With a slight abuse of notations, we will often drop the subscript whenever it is clear from 
the context and write j^-, ■ instead of j^-, ■ , i = 1,2. With these notations, one can write 
the dynamics (jl.ip only in terms of the empirical measure (jl.GI) 

d9i^t = c{9i^t,^^i)dt + ^T^, i9i^t,u}i,Xi)dt+ dBi^t, te[0,r] (2.4) 

where stands for (12.2|] with the choice of g(9,u},x,9,Ld,x) := 

r{9,uj,9,Cj)'df{x,x) and A ;= VN^t- Note also that the McKean-Vlasov equation (jl.7p may 
be written as 

, f) = {J^o, f) + {I's , L[us]f) ds, (2.5) 

where, for any measure A on X x Y x S, the generator K[A] is given by 

^[A]/(t):= ^Ae/(T) + Ve/(r)-|c(6',w) + [r4',A (r)|. 

2.1. Decomposition of the fluctuation process. 


( 2 . 6 ) 


Two-particle fluctuation process. The main difficulty in the analysis of (jl.9p comes from 
the singularity of the spatial kernel 'k (II.3p . In particular, if one studies the process 
alone, one would need to consider test functions {9,uj,x) 1 -^ f{9,L0,x) with singularities 
w.r.t. X whereas the embeddings techniques used in this work require a minimal regularity 
on the test functions (see Section 13.31 below). Hence, the main idea is to introduce an 
auxiliary process TLn va order to bypass the lack of regularity of the kernel 4'. This process, 
that we call two-partiele fluctuation process, is the key object in order to understand the 
influence of the positions {xi,Xj) of the particles on the fluctuations of the whole system. 

Definition 2.1. Define the two-particle fluctuation process by 

:= aiv ( —^ E -i^Ar,t0t't(4'-) 1(^1-, df), te [0,T] (2.7) 

\|A-7v| ijeAjv / 

that is, for any test function {9,uj,x,9,Cb,x) = (r, f) ^ g{T,f) 

, 5) - I 77^ E - 77 ^ E f ^ixi,x)g{Ti^t,f)ut{df) 

\\X^n\ i,jeAN ieAjv 

( 2 . 8 ) 

The process TIn captures the mutual fluctuations of two particles {9i,9j) instead of 
one. One can easily see that TLn captures the correct fluctuations induced by the space 
variables (especially in the supercritical case a> taking g = 1 in (12.81) . one obtains that 

{TLN,t, 1) = which, by LemmaEa 

below and Definition [T3l is exactly of order 1 when a > |. 
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What makes the use oH-Ln critical is that it enables to separate the issue of the sin¬ 
gularity of the spatial kernel 'h from the issue of the regularity of the test functions g: 
it is the process T-Ln itself that carries the singularity in (x,x), not the test functions. 
Hence, we are allowed to (and we will in the sequel) consider test functions as regular as 
required in all variables {9,uj,x), which is crucial for the Sobolev embeddings techniques 
used in the paper. Note that it also necessary to consider the process T-Ln in the subcritical 
case a < I even if ^ = limTv^oo'HAr does not appear in the hnal convergence result (see 
Theorem 12. 7p . 

Relations between r]j^ and RLn- It is immediate to see from (jl.9p and (12. 7h that for all test 
functions (r, f) g{T,f) 

g) = {r]N,t, {’^N,t,^g)), t€[0,T], (2.9) 

where by {r]N,t, ^ 5 )) we mean (??Ar,t(df), {nN,t{dT), ^{x,x)g{T,f))). A natural 

question would be to ask if an equality similar to (12.9p holds in the limit as N ^ 00 

{T-l-t, g) = {r]t, {i't,^g)), te[0,r], (2.10) 

for any possible limits gt and 'HN,t -^n-*oo 'Ht- Supposing that (I2.10h is 

true, a consequence is that the limiting process % becomes a posteriori useless for the 
determination of g: we will show in Theorem 12.71 (see also Remark I4.10p that, using 
(|2.10p . one can characterize the limit oi g^ as the unique solution to a linear SPDE (I2.,'16h 
involving only g and not the auxiliary process Ti. 

Equality (|2.10p is certainly true in the case without space (that is the case considered 
in [ini [27], equivalent to a = 0 in (11.11) 1. One result of the paper is to show that (|2.10l) 
remains true in the subcritical case 0 ^ a < | (although the proof for this equality is not 
straightforward, see Section [T4l and Proposition 14. 9p . 

On the contrary, we will show that (j2.10p does not hold in the supercritical case a > |: 
the limiting processes g and % found in Theorem 12.81 lead to different expressions in 
(|2.10p : when g = 1 m. (|2.10p . one obtains that l(Ht , 1) = {T~Lo , 1) = x{^) ^ 0 (see (|2.38p 
and (|4.3I) 1. On the other hand, {ut, '!'(■, 5)) = /§'k(x,x)dx is a constant C (equal to 
recall Lemma [3lTp and we see from (I2.38P that {gt, C) = 0. 

The reason for this difference is that, when a > |, the predominant scaling (i.e. N^“") 
comes from the rate of convergence of the discrete Riemann sum EjeAjv d'(xi,Xj) to¬ 
wards the integral fg^(xi,x)dx ('Lemma 13.41) . This scaling is strongly related to the 
singularity of the kernel T: if we had replaced 'L by a regular (e.g. C^) function, the 
rate of convergence of this Riemann sum would become and the predominant scaling 
of the whole system would still be \fN. If one naively replaces in (12.91) the empiri¬ 
cal measure by its McKean-Vlasov limit u, the singularity in space disappears: the 
function x ^ [uf, ^{■,x)g{-,f)) is continuous (provided g is). In particular, in (|2.9p . 
the term {gN,t, (ut, 'd/g)) does not contribute to the scaling It only appears in 

{gN,t ) {^N,t - ut, Tgr)), which is a function of 'HN,t but not of gN,t alone. Hence, there is 
no hope to have a closed formula for g when a > |: we obtain as V ^ 00 a system of 
coupled deterministic equations in {g,'H), see (I2.38D . 

Semi-martingale representation of the fluctuation processes. The starting point of the anal¬ 
ysis is to write a semi-martingale decomposition for both processes gjy and T-Ln- We see in 
(j2.11D the use of the process TiN- the singular part in the semi-martingale decomposition 
of r/ 7 v is completely expressed in terms of Rn- 
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Proposition 2.2 (Semi-martingale representation of For every testfunetion {9,u,x) 

f( 6 ,L 0 ,x), for all t e [0,T], one has 

{VN,t, f) = {r]N,o, f) + {r]N,s, L[us]f) ds + (^Ar,^ , ^>[/]) ds + (2.11) 

where L^Vs] is the propagator defined in (j2.6jl and 

= HfW,<^,x,9,uj,x) ■■= S/ 0 f{ 9 ,u;,x)-T{ 9 ,uj, 9 ,oj). 
and where the martingale term is given by 

E te[0,T]. 

’ i^A.- -^0 


( 2 . 12 ) 


(2.13) 


Proposition 2.3 (Semi-martingale representation for T-Ln)- For any regular and bounded 
test function (r, f) ^ g{T,f) 

, g) = {Td-Nfi , g) + {'Hn,s , ^sg) ds + i: FN,sg ds + i: GN,sgds + g, (2.14) 


where 

for 


^sg{T, f) ;= f) + ^G)g{T, f) 


■^y^g{'r,f) := -Agg 5 (r,f) + Veg{r,f) ■ {c{9,uj) + [r4',i^s] (r)} 

+ V 0 fif(r,f) ■ {c{9,uj) + (f)} , 

■^s^^g{x,f) ■■= {us, ^{■,x)V0i-,T)) ■T{9,uj,9,oj), 
where the remaining terms F^ and Gn are 

FN,sg ■= {'Hn,s , V05(r,f) ■ [r4',j/Ar,s - i^s] (e)) 

GN,sg ■= {'Hn,s , {t^N,s - I's, 4'(',x)V05'(-,t)) ■ T{9,u,9,u)) 


(2.15) 


(2.16) 

(2.17) 

(2.18) 
(2.19) 


('H) 

and the martingale part A4)y f is given by 

E f ^§9iF,s,Tj,s)^{Xi,Xj) 

A^r ocA .. -'U 


l^ivl ieA, 


IAa^I ijeAjv 
't(j^ 




E f ^77^ E ^egin,s,Tj^s)'^{xi,Xj) - f V9g{Ti^s,f)'it{xi,x)us{df) 

= A Jo \ Aat .. j I 


dB 


l^S 

( 2 . 20 ) 


The proofs of Propositions 12.21 and 12.31 are given in Section 15.11 The whole point of the 
paper is to take the limit as A ^ oo in the semi-martingale decompositions (I2.11|) and 
(TOD . 

2.2. Assumptions. Define the following integer 

P := m + n + 1 = dim(X x Y x S). (2-21) 

The local dynamics term {9,uj) c{9,lj) in (jl.lh is supposed to be differentiable w.r.t. 

{9,u}) up to order 3P + 9 and satisfies a one-sided Lipschitz condition w.r.t. 9, uniformly 

sup{(0-0)-(c(0,w)-c(0,a;))} ^L\9-9f, (0,0) eX^ 

ojeY 


m (w: 


( 2 . 22 ) 
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for some constant L. We also suppose that i-s- c{9,lj) and its derivatives up to order 

3P + 9 are bounded in 9 with polynomial bound in uj: there exist t ^ 1 and C > 0 such 
that for all differential operator D in {9,u}) of order smaller than 3P + 9 

sup|L)c(0,u;)| ^C(l + |u;|''), w e Y. (2.23) 

6»eX 


The interaction term T in (jl.ip is supposed to be bounded by || T ||^, globally Lipschitz- 
continuous on (X x Y)^, with a Lipschitz constant || T 11^^^ and differentiable in all variables 
up to order 3P + 9 with bounded derivatives. We also require that 


E sup( [ \Dl^T{9M0,u;)\d9) 
|fc| ' 

Define the following exponents (where l is given by (12.231) 1 


< +0O. 




7 := max yP, 

K:=m + 1 and k'-= k + 2'y, 
t := n + 2i + 1 and T ■= l + 2'y. 


(2.24) 


(2.25) 

(2.26) 
(2.27) 


The law ( of the particles (0i,o)ieAjv is assumed to be absolutely continuous with respect 
to the Lebesgue measure on X = and its density (that we also denote by 9 C(^) with 

a slight abuse of notation) satisfies the integrability condition 

3p>m, ^C(0)^cl0<+oo. (2.28) 

We also assume that C and satisfy the moment conditions 

y' C( d0) <+ 0 O and /r( dw) <+oo, (2.29) 

where (R, t) are given in (j2.26p and (|2.27|) . In what follows, we denote as £(■) the expecta¬ 
tion with respect to the initial condition 9ifi, the Brownian motions {Bi) and the disorder 
(cji) (i.e. we work in the averaged model). 


2.3. Estimates on the McKean-Vlasov eqnation. We first recall some results con¬ 
cerning the continuous limit (|1.7p . The crucial object here is the nonlinear process as¬ 
sociated to the McKean-Vlasov equation (fTTHl [Ml [28]. The following result is a direct 
consequence of [28]) Section 3.2: 


Proposition 2.4. Under the assumptions made in Section [2.2l there is path'wise existence 
and uniqueness of the solution (^^’^,a;,x)jg[o,T] of the nonlinear system 

= + £ (c{9f’‘^,uj) + [tT, A,](0f’",a;,x)) ds + Bt, 

Xt = X^’‘^{d9)fi{duj) dx, te[0,T]. (2.30) 

X^’‘^{d9)p{dLo) is the law of {9^’‘^ ,lo), 

satisfying sup^^-p / v |a;|®'') Ai(d0, dw, dx) < +oo, where R and I are given in (I2.26h 
and (I2.27p . 


Proposition 2.5. Under the hypothesis made in Section \2.2l there exists a unique weak 
solution t^ ft in C([0,r], A4i(R”^ x R”)) to (|1.7p satisfying 


sup 

t^T 



V nt(d9, duj, dx) 


< + 00 . 
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Moreover, there exists a continuous measure-valued process t on X x Y such that 


t't( d0, dto, dx) = d0, doj) dx, te[0,T]. (2-31) 

The process ^ admits a regular density {t,6,uj) pt{9,u) with respect to d0 0 /i(da;) 

^t{d.6, duj) = pt{0,u))d0pL{(luj), te[0,T], (2.32) 

and this density p satisfies the a priori estimates: 

0^pt(6*,a;)^ , weY, te(0,T] (2.33) 

1 “H I I ^ 

|div 0 pt( 6 ',w)| ^ weY, ie( 0 ,r], (2.34) 

i“0+2 


for some 0 < oq < | ■ 

The proof of Proposition 12.51 is postponed to Section 15.21 

Proposition 2.6. The particle system (HH) and the non-linear process (I2.30p satisfy the 
moment conditions 


sup sup E I sup \6i^. 

N ^ 1 ieAjsj \s 


\8k, 


< + 0 O and E I sup sup \d' 

\s ^ T xeS 




< + 00 , 


(2.35) 


where k is given in (I2.26h . 

Proof of Proposition \2.6[ The estimate on the nonlinear process 0 is a direct consequence 
of Proposition 12.41 The same estimate for the particle system is standard and left to the 
reader. □ 


2.4. Fluctuations results. We prove the convergence of the fluctuations processes in 
appropriate Sobolev weighted spaces of distributions V and W that are defined in the 
next section, using Hilbertian techniques developed by Fernandez and Meleard in |16j . 
We only state the result here and refer to Section 13.31 for precise definitions of these 
spaces. 

Fluctuation result in the suhcritical case a < |. The first main result of the paper is the 
following (recall the definition of the propagator L[zzs] in (I2.6p l: 

Theoreui 2.7. Suppose a< Under the hypothesis of Section \2.iA the random process 
{vn)n^i defined by (11.91) converges in law as N ->■ oo to rj € C([0,T], ) solution 

in linear stochastic partial differential equation 

r]t = 'no +C*s'nsds + M\'^\ te[0,r] (2.36) 

where 

Csf{0,uj,x) ;= L[vs]f{e,u),x) + {vs , V 0 /(-) ■ r(-,-, 6 ',w)T(-,x)), (2.37) 

go is a Gaussian process with explicit covariance given in Proposition \f. 1\ and is an 

explicit martingale given in Definition \4-2\ go o,f^d are independent. 
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Convergence result when a > |. The second main result concerns the fluctuations of m 
in the supercritical case. For technical reasons (see Section IT3]1 . in addition to the assump¬ 
tions of Section EH we restrict here to the case where X is no longer but a compact 
domain of M™. The example we have particularly in mind here is the Kuramoto case where 
X = M/27rZ. We also suppose that the support of the distribution // of the disorder is com¬ 
pact. These further assumptions are made in order to ensure uniqueness of a solution to 
(|2..S8p below, due to the nonstandard nature of the operator (see Section [45]) . 


Theorem 2.8. Suppose a> Under the hypothesis of Sections \2.S\ and \2.4\ the random 
process {riN,'HN)N ^ i converges in law as N oo to solution in 

^%(p+ 2 )) system of coupled deterministic equations 


'nt = L[ns]*'nsds + ^*'Hsds, 

'Ht = 'Ho+ r ^fUsds, 

Jo 


t ^ [0)^]) 


(2.38) 


where T-Lq is a nontrivial initial condition defined in Proposition \4.1\ and LYvg]* (respec¬ 
tively <1>’^ and .if/ ) is the dual of the propagator LYvs] defined in (j2.6p (respectively of the 
linear form <h defined in (j2.12l) and the linear operator defined in (|2.15l ) ). 


Remark 2.9. The assumptions of Theorems \2.'i/\ and \2.8\ do not cover the case of FitzHugh- 
Nagumo oscillators (that is when only a one-sided Lipchitz continuity on c as in (|2.22l) 
and polynomial bound on c is required). A careful reading of the following shows that the 
tightness results fTheorems \3.23\ and \3.24\) are indeed true in the FitzHugh-Nagumo case. 
The restrictive conditions of the paper are only required for the uniqueness of the limits. 

Comments on the critical fluctuations a = |. In the critical case a = |, it is expected that 
the correct scaling is otv = \/lV although the techniques used in this work do not seem 
to make this intuition rigorous. Nevertheless, a closer look at the proofs below shows the 
following partial result 


Proposition 2.10 (Critical fluctuations). When a = ^, under the assumptions of Sec- 
tion \2.2[ the following convergence holds ^ C([0,T], 


^/N 

iEn 


{uN - p) 


0 , 


as N ->■ oo. 


(2.39) 


3. Tightness results 

This section is devoted to the tightness of t/at and LLn, based on their semi-martingale 
decomposition (|2.11|1 and (I2.14p . 

3.1. First estimates on T. 

Lemma 3.1. There exists a constant C that only depends on a such for all {x,y,z) 

(3.1) 

Proof of Lemma \3.1\ Straightforward (see for example m, Lemma 2.5 for a proof). □ 






















12 


ERIC LUgONi AND WILHELM STANNAT^-^ 


Lemma 3.2. For all /3 > 0, there exists a constant C > 0 (that only depends on (3), such 
that for all N ^ 1, for all i e Ajv, 




< C- 


N i/0<^<l, 

iVlniV if (3 = 1, 
iV^ if^>l. 


j*i 

Proof of Lemm,a, \S.^ This lemma has been proven in [28], Lemma 6.1. 


(3.2) 

□ 


Remark 3.3. An easy consequence of Lemma \3.f^ that will be continuously used in the 
following is that, since a € [0,1), 


sup sup ——■ j: >i-(x.,x,)<+oc. 

N ^ 1 ieAjv jeAjv 


(3.3) 


The following (very simple) lemma is at the core of the difficulties of the paper: the 
rate of convergence, as oo, of the Riemann sum associated to the function x ^ 
to its integral is which is in particular smaller than the Gaussian scaling \/N, when 

a > 2- 

Lemma 3.4. For all a e [0,1), ^ 1 and i e A^, the quantity 

E '^{xi,Xj)-£^{xi,x)dx^ (3.4) 

is independent of i € An and converges as N ^ oo to a constant x(a) ^ 0 (that depends 
only on a). 


Proof of Lemma \3.4\ Prom the definition of d{-,-) in (|1.2I) . a direct calculation shows that 
/g'I'(xi, x) dx is actually independent of z e and equal to Moreover, one directly 

sees from (fL2]l that EjeAjv '^ixi,Xj) = (E^i By a usual comparison 

with integrals, 


N I 


^ 1,1 1 

— dt H-1- 

2 2N^ 




a 


N 


l-a 


1 


1 


1 ^ 

l-a ~ l-a 2 2A1" 


-rnm 


{u + A:)“+1 
N-l 


du, 


a 


Defining C'(a) := | “ l) (n+fc)^+i dw * 0, one obtains 


du. 


N I ^ pjl-a 

hi 1 -« 


+ C{a) + 


2m 




a 


k^Niu + k) 


OL+l 


du. 


Consequently, 


1 2 
^ E '^{Xi,Xj) = — 


ATl-a 

2 “ 


/Ari-« 

\ l-a 


+ C{a) + rN 


) 


2 “C(a) 2" 

+ —rN - 


20-1 

20-1 


l-a N ’ 

where ttv 0 as Ai -> oo. Lemma 13.41 follows with x(a) := 2"C'(a). 


□ 
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3.2. Estimates on the nonlinear process. The tightness result is based on a coupling 
argument: for all i e A^r, define the nonlinear process 6i associated to the diffusion 9i, 
with the same Brownian motion Bi, initial condition Oi^, disorder oji and position Xji 


di,t = 0i,o + {6i^s,uJi,Xi)^ ds + Bi^t, te[0,r], zeAw. (3.5) 


Remark 3.5. Note that, due to the rotational invariance of§, the function {6,uj,x) 
[TTjt't] (9,uj,x) is actually independent of x (see (15.81) ). so that the definition of 9i (con¬ 
trary to the the microscopic particle system 9i^N in (|2.4p ) does not depend on the position 
Xi- But we keep this dependence for consistency of notations. 


The following proposition quantifies how good the approximation of 9i by 9i is: 


Proposition 3.6. Recall the definition of the scaling parameter a n in (jl.lOp . There exists 
a constant C > 0 only depending on T, c, T and T such that 

E [max sup ^N ^ 1. (3.6) 

\ieANt^T ) 

Proof of Proposition 13.61 is given in Section 15.31 


3.3. Weighted-Sobolev spaces. We introduce weighted spaces of distributions [161 [27] 
for the study of processes ryw and TLn- The definitions are given in two steps, since ryw 
acts on test functions of one variable r ns- /(r) whereas TiN acts on test functions of two 
variables (r, f) g{T,f). The only thing that differs from the usual Sobolev norm is 
the presence of a polynomial weight made necessary for the control in (jl.ip of the term 
{9,lj) i-s- c{9,ijj) on X X Y. 


Sobolev spaces on X x Y x S. Recall the notation r = (9,uj,x) e XxY x§ and the definition 
of K, i in (|2.23p and of P in (|2.2ip . For all 1,1' 0 and pfiO, define the norm || ■ ||p ^ ^ over 

test functions r f(T) by 


/ 




s / 




■ dr 


(3.7) 


JA:| (1 + fit + / 

where, if k = {ki,... ,kp) and r = {9,io,x) = {9^^\... ... ,io^'''\x), we define 

\k\ ■■= T,fLi h and D'ff{T) := Denote by the com¬ 
pletion of regular functions with compact support under the norm || ■ ^, || ■ ^ ^ j 


/ I' 

is an Hilbert space. Denote by Vjp^ its dual, with norm || ■ ||_p; p. 

I I' 

Denote as Cf the Banach space of functions / with continuous derivatives up to order 
p such that lim| 0 |+|t^|^oo sup^-^g = 0, for all \k\ ^ p. Endow this space with the 

norm || • given by 


IIC*’ 


f\\c‘^i'-= E sup 


\Drf{r)\ 

^ p 7 1^1 |u^l 


One has the continuous embedding (see [ 2 ] or [16] Section 2 . 1 ): 


V 


1 , 1 ' 


q+r 


&/, rfid, q>Pl2, 1,1'fiO, 

^^hh' , Pn , ,, n ,, 

Yfi , r 0, k > — +l k > — + l , 
2 2 


(3.8) 

(3.9) 
(3.10) 
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and the Hilbert-Schmidt embedding 


v: 


1,1' 

q+r 


yUk,l'+k^ r ^ 0, q> P/2, k > P/2, I, l' ^ 0. 


The corresponding Hilbert-Schmidt embedding holds for the dual spaces: 

rl+k,l'+k 


v_ 


V' 


-(q+r 


y r ^ 0, g > P/2, /c > P/2, l,l'^ 


0 . 


(3.11) 

(3.12) 


Sobolev spaces on (X x Y x S)^. We now define similar Sobolev spaces for test functions of 
two variables (r, f) g(T,f). For all 1,1' ^0 and p ^ 0, define: 


/ 


I 9 \\p,i,i' 


S / 

|fc|+|/c| ^ p 




1 

\ 2 


and 


9 Ll*' - E sup 


w{T,f,l,l'Y 

P^P|(7(T,f) 


dr dr 


/ 


|fc|+|fc| 




where the weight w is given by 

w{T,f,l,l') := ^1 + \9f + ^1 + |a;|^ + \Cj'( j . 


(3.13) 


(3.14) 


(3.15) 


Define similarly || ■ ||p ^ j, (wi:E II ' ILp,;,/') and {Cp’' ,\\-\\^i,i') the corresponding 

spaces. We use for simplicity the same notations for the norms in Sections 13.31 and 13.31 
since the distinction between both cases will be often clear from the context. Similarly, 
one has the embeddings: 


WqE cl'’' , r ^ 0, q> P, l,l' ^ 0, (3.16) 

&/ -+ r ^ 0, k>m + l, k'>n + l', (3.17) 

^q+r ^ r ^ 0, q>P, k>P, 1,1' ^0, (3.18) 

^l_+k,l'+k ^ wC^^^^y r^O, q>P, k>P 1,1'^ 0. (3.19) 


For the proof of tightness of (i^atjPiv) and the identification of its limit, we will need in 
Section 0] several instances of the spaces V and W, for different choices of the Sobolev 
parameters {p,l,l'). In the rest of Section [3l we work with general parameters and the 
precise values of the parameters will be specified in Section 01 


Definition 3.7. Fix two integers kq ^ 0 and lq ^ 0 and define (recall (I2.25P ) 

(fi:i,ii) := (kq+ 7Po+7)- (3-20) 

Note that the definition of (kq,ki, ioPi) ensures that the Sobolev embeddings (13.lip . 
P3.12p . p3.18n and (I3.19P are true with I = kq, I' = lq and fe = 7 (provided the regularity 
index q is such that q> P). 


Continuity of linear forms. We place ourselves in the context of test functions of two 
variables (Section 13.3p . Define the following linear forms describing the variations of the 
test functions with respect to either {9, to) or to the space variable x. 
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Definition 3.8 (Continuity in 9). For fixed r = {9,uj,x), f = (0,u),x) e X x Y x and 
7r,7r e X, define the linear forms: 

'R,r,r, 7 T,n{g) ■= g{9,u!,x,9,oj,x) - g{Tr,uj,x,TT,uj,x), (3.21) 

<ST,r{g) ■■= g{T,f), (3.22) 

TT,f{g) ■■= div0g{T, f) + divg5f(r, f), (3.23) 


Definition 3.9 (Continuity in x). For any a,b € E>, A subinterval of§, r e X x Y x S, 
N I, i € An, (n,... ,tn) e (X X Y X E>)^, r = I,... ,m, ^ e Ali(X x Y), define the linear 
forms 


U{g) = Ur,a,b, A,da) ■= f {9{r,0,Cj,a) - g{T,e,dj,x))C{de, da;)dx, 

Vig) = VN,i,Tig) ■■= 7T^ E dg(r)g{Tk,Ti)^{xk,Xi), 

I^Af| fceAjv 

W{g) = WN,i,r,d9) ■= J dg(r) gin, f)^{xi,x)C{de, du))dx. 


(3.24) 

(3.25) 

(3.26) 


Proposition 3.10. For all r, f,7r,7r, for any P + 2, the linear forms of Definition \3.t^ 
are continuous on with norms: 

II 7^r,f,7r,s■ Lg,nuii ^ C* (1^ - 7r| + |0 - 7r|) (r, f, vr, vr), 

ll>5r,f ILq,Ki,ti ^Cw{T,f,Kl,tl), 
where w is defined in (I3.15P and 

XKi,u('r,'?,7r,7r) := 1 + l^l'^^ + l^l^^i + |7r|''^ + |7r|'"i + |a;|'^ + (3.27) 

and where the constant C is independent o/(r, f, vr, vr). 


Proof of Provosition \3.1(K We prove the first estimate of Proposition 13.101 Fix g regular, 
with compact support in {9,uj,6,uj). 

m 

|^r,f, 7 r,s-( 5 ')| ^ (|6»-0| + |7r-7r|)sup Y,{dr,(r)g{u,u),x,v,Cj,x) + d^(r)g{u,uj,x,v,Co,x)) , 

U^V 

^ C (l^ - 0| + IvT - 7f|) Xkuh (a t", vr, vr) || g ||^-i,h , 

^ C (10 - 0| + |vr - t|) {r, f, t, it) || g , (by UJM) 

where the supremum in the first inequality is taken over |u| ^ |0| + |0| and |u| ^ |vr| + |vr| 
and where the integer P is defined in (|2.21l) . The result follows by a density argument. 
The proofs of the estimates on S and T are similar and left to the reader. □ 


Proposition 3.11. For any P+ 2, the linear forms of Definition \3.!^ are continuous 
on in the sense that 


l^ll-qAi,.i ^ C(l + |0|''i + |a;|‘')^sup|a-x|J^T(6,x)dx, 

c 

II ll-l},Kl,tl ^ I A I E w{Tk,Ti,Kl,il)^{Xk,Xi), 

fceAjv 

'^\\-q,Ki,n ^ ^ J w{f,Ti,Ki,Li)f,{de, dw). 


(3.28) 

(3.29) 

(3.30) 
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Proof of Proposition I,?. 1 A Using the moment estimates of Proposition 12.61 
P{g)\^ j^^{b,x)^j{g{T,e,Cb,a)-g{T,e,u},x))i{de,(lu}) dx 
^ C II 5 II 'ID'Tj ^ dcxj) ^X) I fl X I dx; 


^C\\g\ 


q,K,i,Li 


( 


(1 + l^l''^ + |a;|''^) sup|x - a| 1 / ^'( 6 ,x)dx. 


^,xeA / J^ 

The proof for V and W are similar and left to the reader. 


□ 


3.4. Tightness criterion. We prove tightness for both fluctuation processes r/jv (ll.9|l and 

P.N pup in suitable V and W defined in Section 13.31 We recall the tightness criterion 
used here ([ 22 ], p.35): a sequence of (OAr,.TAr_t)-adapted processes with paths 

in C{[0,T], H) where H is an Hilbert space is tight if both conditions hold: 

(1) There exists an Hilbert space PIq such that Hq ^ H (with Hilbert-Schmidt injec¬ 
tion) and such that for all t ^ T, 

supE(|| Yat,! 11^ ) <+ 0 O, (3.31) 

N “ 

(2) Aldous condition: for every £ 1,62 > 0, there exists sq > 0 and an integer Nq such 
that for every (JT/v_t)-stopping time Tn ^T, 

sup sup P(|| Yn,tn - YN,rN+s \\h ^ ^i) ^ £ 2 - (3.32) 

N ^ s ^ so 

3.5. Tightness of the two-particle fluctnation process. 


Boundedness of the fluetuation proeess. 

Proposition 3.12. Under the assumptions made in Section WfE. for any q ^ P + 2, 

sup sup E (II HN,t \\t„ Ki i, ) < (3-33) 

Remark 3.13. We deduce from (j3.33p f/iaf sup;^ ^ < t^ (ll ll-q ki ti) which 

is the main estimate that we use in the following. A fourth moment estimate is only needed 
in Provosition \3.1(A 

Proof of Proposition 13.121 is given in Section 15.41 


Proposition 3.14. Under the assumptions of Section \2.Sl the linear operator defined 
in ()2.15l) is continuous from to there exists some constant C > 0, such 

that for all g € 

\\-^sg\\q^Kl,Ll ^ ^ II S' llg+P+2,K0,t0 • (3.34) 


Proof of Proposition \3A^ The first term ^Ag^g in A^'gg clearly satisfies (13.341) . We con¬ 
centrate on the second term (r, f) Veg{T,f) ■ c(9,lj): since || Vgg{T,f)c{6,u!) ||g^^ = 

\D!;DlWeg(r,fM0,uj))f 


■'j/cl + lfcl : 


,1 


— 7 —=- ^ -^drdf, it suffices to estimate the quantity 

|^iV05f(T^)l^l-D2c(6',a;)|" 


/ 


w{T,f,Ki,Liy 


■ dr df, 
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for all differential operators Di and D 2 of order smaller or equal than q. Namely, 


/ 


\Diye9{r,f)\^ \D2 c{0 

w{T,f,Ki,iiy 


dr df 


|^iVec/(r,f)| \D 2 c{e,u))\ w{T,f,KQ,iQf 


w{T,f,Ko,LQy 


-! 

^ II II 2 r 

^ 9+1 J 




dr df, 


\D2c{e, uj)f' w{t, f, Kq, io)" 


u;(r,f,Ni,ii)2 


dr df, 

(3.35) 


II II 2 

^ ^ II 5 llg+P+2,K0, 


^•0 


thanks to the embedding ()3.16p and by assumption on c (recall (12.221) 1. Note that the 
definition of 7 in ()2.25p ensures that the integral in (I3.35P is indeed finite. The third term 
in p2.15p . (r, f) Ve 5 '(r, f) ■ [rTji/s] (r) can be treated in the exact same way, observing 
that [T'hjUs] (t) is actually independent of the space variable x (see (j5.8p l. It remains to 
treat the last term (r,f) j 0 g{-,T)) ■ T{9,uj,9,oj). For all r = {9,uj,x), 


{us,'i>{-,x)V 0 g{-,T)) = /^(y,x) V 0 g{TT, 4 >,y, 9 ,u},x)^s{dTT, d(j))dy, 

= P{z) / Vng{TT,(j),z + x, 9 ,u},x)^s{d 7 r, d(j))dz, 

w/§ -^XxY 


where (recall (II.2p and (ll.3jl ) 


p{z) ■■= min(|z|, 1 - | 2 ;|) “ 


(3.36) 


belongs to L^(S). In particular, for any differential operator D (acting on {9,u!,x)) of 
order smaller or equal than q, 

\D{us,d>{-,x)S/Qg{-,T))\= [p{z) f D [v gg{'n:,cl),z + x, 9, uj,x)]^s{ dir, del)) dz, 

JE> JXxY 

< c II 5 iic"0’‘o II p IILi(s) f (1+Kr°+i^r°) ( 1 +i<^r°+i^^r®)6(d^, 

<c'ii9ii,.p,2,„...(i+i9r + M‘»). 

since has finite moments of order kq in 9 and lq in oj (recall Proposition 12.61) . Since 
r(',-) is supposed to be regular and bounded as well as its derivatives (recall Section 
it suffices to estimate, 


/ 


D{i 2 s, T(-,x)Vgg(-,r))f 

w{T,f,Kl,biy 


dr df ^ C 


g 


2 

g+P+ 2 ,Ko,to 


/ 


(l + | 6»|'‘0 + | w |‘ 0 )^ 

u;(r,f,Ki,ii )2 


dr df, 


<C\\g 


II 2 

IIij+P+2,ko,'-0 ■ 


This concludes the proof of Proposition 13.141 


□ 


Proposition 3.15. Under the assumptions of SectionlKM for any q ^ P + 2, <h defined in 
P 2 . 12 p is a continuous linear operator from ^ql’p ^2 i.e. there exists a constant 

C > 0 such that for all f e 

II d’[/] llq,Ki,ti ^ C' II / ||q+p+2,Ko+to ■ (3.37) 

Proof of Proposition \d.l5l Straightforward, since by assumption, P is bounded as well as 
its derivatives. □ 
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Proposition 3.16. For any q ^ P + 2, the remaining terms Fjsf and Gm in the semi¬ 
martingale decomposition (I2.14|) of FLm converge uniformly in time to Q, as N ^ oo, in 


more precisely, there exists a constant C > 0 such that 


sup E II Fn, 


t<T 


\ c 

H—> 


-q,Ki,Li 


supE(|| Gjv.i ^ 

t sg T ' i. 1/ 


t s£ T 


(3.38) 

(3.39) 


Proof of Proposition Writing explicitly Etv and Gjv from their definitions (I2.18[) and 

(j2.19p gives: 

FN,tg = 7T^ E nN,t - nt] (Ti^t) 

I^A^I kAn 

(3.40) 


aiv 


E ^{xi,Xj)Vegipt,Tj,t) - f '^ixi,x)V0giTi^t,f)i't{df)]\, 

I^AfljeAjv I) 


GN,tg=T^^ E {^N,t-nt,'^{-,Xi)VQg{-,Ti^t)) 


(3.41) 


For all r = (0,a;,x) e X X Y X §, define the linear form pN,T,t by 
\^N\j,A, 


^N,T,tg •= 0 ,N 

Hence, 


E ^{x,Xj)Veg{x,Tj^t)- j ^{x,x)Veg{T,f)vt{df)Y (3.42) 


E E(|[r4/,i.^,i-i.i](r,,t)h^E(||^^,.,,,El' )"> 

I^A^I z^An ^ 

by Cauchy-Schwartz inequality. By the same procedure as in Proposition 13.121 in Sec¬ 
tion EH (note that the structure of pN,Ti t,t is very similar to Pn, recall (12.81) 1. it is easy 
to show that 

c 

sup sup E(|[r^,i/Ar,t - i/j] (ri,t)|'^) ^ , 

t zgAtv 

II l|4 

sup sup sup E( 

A 3s 11 sg T ieAiv ^ 

Hence, ()3.38D easily follows. The same argument holds for Gn. Proposition 13.16) is proven. 

□ 

/'^'V 

Proposition 3.17. For any q ^ P + 2, the process Al)^ ' defined in (I2.20p is a martingale 
with values in and there exists a constant G > 0 such that 


sup E I 

t<T 




I^aI 


(3.43) 
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Remark 3.18. Due to the embedding (I3.19p . one has || • ||_(q+p+ 2 ) kq to ^ ^ II ’ H-g ki n’ 
that every estimate (especially (|3.38|) . (|3.39p and (|3.43p ) involving the || ■ Il_g 
also valid for the || ■ ||_(q+p+ 2 ) kq 


Proof of Proposition \3Al\ For any test function /, 






ON 


I-^nI iJeAff 
UN 


rt 

E / ^§9{Ds,Tj^s)'^{Xi,Xj)- 

icA “'O 


dB 


3,s 


l-^Afl 


E f yeg{Ti^s,Tj,s)^{xi,Xj) ■ dBi^s 
lA I E f f V8g{Ti^s,f)'^{xi,x)us{df) ■ 

I^A^I teAjv 


ON 


dB 


2,S 




( 2 ), 


<(3). 


We only treat the first term the other two being similar: using a complete orthonor¬ 
mal system (</?p)p ^ i in (we drop the dependence in ki, below for simplicity), we 

have successively 


E sup 

\t^T 


M 


m 

N,t 


]=e(sup S < E e( 

V 1 j p^ 1 V 


sup 

tCT 


-^Ar,i Pp 


1 




EEf ■^ 7 \rr^p| ] ’ Doob’s inequality), 

p > 1 ' ’ ' 






IAat] i,j,keArf 
'T 

\P 


Y, ^{Xi,Xj)'^{Xk,Xj) 


e( ^ |Vg95p(Ti^s, Tj^s)! j e( ^ Q^Pp{Tk^s,Tj,^)\ I 




‘‘N 


I-^nI i,j,keA 


^ 'if{xi,Xj)^{xk,Xj) 




(3.44) 




/~T 

^77^ E '^{Xi,Xj)^{Xk,Xj) / 'E{w{Ti^s,Tj,s))'E,{w{Tk,s,rj^s 

IAatI i,j,k€AN 

(3.45) 

(3.46) 






E '^ixi,Xj)^{xk,Xj) 

|Ajv| i,j,fceAjv I^AtI 


)) ds, 


where we recall in (I3.44p the definition ^3.231) of the linear form T, where we used Propo¬ 
sition 13.101 in p3.45p and Lemma 13.21 in p3.46p . This concludes the proof of Proposi¬ 
tion 13.171 □ 
























































20 


ERIC LUgONi AND WILHELM STANNAT^-^ 


Proposition 3.19. Under the hypotheses of Section \2. 21 for any q ^ P+2, the two-particle 
fluctuation process T-Ln belongs uniformly to 

i^up E (sup II 'HN,t ll!(g+P+2),^o/o) < 

Proof of Proposition \2.1fA Let {'i/jp)p^i be a complete orthonormal system in W^°p‘|^ 2 - 
Then, for all p ^ 1, 


, ^Ppf ^ cl {TiNfi , f^pf + t£ {Pn,s , -^sfjp? ds + T \Fn,s'^p? ds 


+ T 


IGat^si/’pI^ ds + M^^Jifp 


so that, by Doob’s inequality. 


e( ^ sup(?{7v,t,i/>p)M ^qEl ^ {nNfl,i^p?\+T r^El ^ 

\p^it^T j y \p>i / \p>i 


,f]ds 


siApP ds 


+ E ^ 

\p^ 1 


Since E (sup^ ^ 7 . || ||^(g^p_^ 2 ),Ko,to) ^ E (llp ^ ^ sup^ ^ , lAp)^), we obtain from 

the previous bound that 




)<c(e(| 


2 

T~Ln,0 ll-(g+P+2) 




X’p)' 


+ T 


+ EI 


) ds 
3 / 


E(||FAr,s||_(,+p+2),Ko,to, 

' )]• 

-{q+P+2),Ko,i.o / I 


M 


(H) 

■N,T 


(3.48) 


Using ()3.43l) and Remark 13.181 one obtains that 


sup E I 

1 ^ A V 


M 


(n) 

N,T 


3:''o) 


< + 00 , 


-(g+P+2),Ko,, 

and using (I3.33p . one obtains a similar bound for the initial fluctuations 

^sup e(||?^^,o||Vp,2),.o/o) 


< +00. 


(3.49) 


(3.50) 
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Moreover, if one introduces the linear form {T~{-n,s , •^s'4^)^ we see from Proposi¬ 

tion [HTH] that iN,s is continuous on W(q+p+2),reo,to- iiideed, for all ij: e 'W(q+p+2),K,o,i-o 

\^N,s'>P\ < II 'Hn,S ILg,Kl,tl II llg,Kl ,0 ’ 

< c II 'Hn,s II i’ Il 5 +P+ 2 ,K 0,^0 ’ (thanks to ^M)- 

In particular 

E ( (Wjv,. . = E (I II?,, 

<Ce(||H«p||?,„„), 

^Csup supE(||'HAr, 5 ||^ )<+oo, (by (|3.33|)). 

is:iVss:r ' i, 1/ 

Putting the previous estimate as well as (j3.38D . (I3.39|) . (|3.49p and (|3.5Up into (I3.48D . one 
obtains the result (|3.47p . Proposition 13.191 is proven. □ 

The rest of this section is devoted to prove the tightness of the two-particle fluctuation 
process Hn- We follow here the lines of HU [231127!. 

Proposition 3.20. For all q ^ P + 2, for every N 1, the trajectories of the two-particle 
fluctuation process Pn o,nd the martingale are almost-surely continuous in W^“p'(^ 2 - 

Proof of Proposition WjWK We only prove the result for Pn, the proof for ^ being 
similar. Using again a complete orthonormal system {'f’p)p^ i in we know that, 

thanks to the proof of Proposition 13.191 with probability 1, for all e > 0, there exists po ^ 1 
such that Ep>po sapt ^ t , fl'p)^ < Let 0 ^t and tm a sequence converging to 
t, as m ^ oo. Then, 

II - '^N,t ll-(g+p+2),«;o,to " S - '^N,t , f^pf , 

p ^ 1 

Po 

^ S , f^pf + 2 XI {{'^N,trr, ) V’p)^ + , f^pf) , 

p-1 p>po 

PQ 2e 

^ X - ^N,t , Ippf + — ^ e, 

p=l a 

for tm sufficiently close to t (by continuity of the map 1 {PN,t , tpp) for all p ^ po). □ 
Proposition 3.21. The process P^ is a semimartingale in this space 

PN,t = Pn,o + J'^ ■^s'^N,s(^S + J'^ FiSf^sds + i: GN,sds + M^^l (3.51) 

/'jj\ 

where is the adjoint of and Ai)^ is a martingale in Doob-Meyer 

process ((AI^^)) with values in is given for every tp, ip ^ ^q+’p +2 

by 

((AI!?^))i(7’)(V') = bN,s{’^N,s){<f,'fl^)ds, (3.52) 
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where for any measure A onXxYxS (we recall (12.3|] and write 9 = {9^^\ ..., e XJ, 


m p 

£ / [(9g»g))^',A]2(f)[(<9e»iA)^',A]2(f)A(df) 

r=l 

m r ( \ 

r=l 

m r / N 

r=l 

r=l~' ' 


and 


(3.53) 


•= E - / (9e(r)g>(r,f)^'(x,x)i/s(df). (3.54) 


Proof of Provosition [071 It is a simple adaptation to our case of [16], Proposition 4.5. 
We refer to it for details. Note that, using the independence of the Brownian motions 
one has 


= 77^ E E f dg(r)p{Ti^s,Tj,s)'^{Xi,Xj)dg^r)fj{Tk,s,Tj,s)'^{Xk,Xj)ds 

* IAatI ij,A:eAjvr-=l-^0 




m . s 

EE. de(.-)P{p,s, Tj,s)'^{xi,Xj)e)(’^{Tj^s, ds 

;eA«r=l'^0 


IAat] ijeAjvr'= 


A AT 


™ r* 1) 

EE. Tj,s)4'(xi, Xj)e%'^{Tj^s, p) ds 

i,nh-N r=l 


lAivr 


i,jeAN 


IAatI^ ifX 


E E £ ^Nl(p,s,p)eP{Ti^s:'ilj)ds, 

EA„r=l 


which is precisely (|3.52l) . 


□ 


Tightness results. 

Proposition 3.22. The sequence of the laws o/(A4^^)7v is tight in c([o,t],w::°^^;^^2))- 

/ TJT\ 

Proof of Proposition \d.2‘A We prove that ^ satisfies the tightness criterion of Sec¬ 
tion 13.41 by Proposition 13.171 the first part of the tightness criterion above is verified for 
Hq - Wlq’‘'^ and H - > since the embedding ^ Hilbert- 

Schmidt. The second part of the criterion is verified if (by Rebolledo’s Theorem, see |22] . 
p.40) it is satisfied by the trace tr^^o.^o ((AIl^ A). Let ei > 0 and Tjv ^ T be a stopping 

-(.q+P+2) 
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time. For a complete orthonormal system {'ipp)p^ i in have 

P ( ((A^^ ^)}tjv+s ~ ((-^AT ^)}Tiv ^ ^i) ) 

\ -(<J+e+2) -(g+P+2) / 




-Hr 

£l \JTi^ 


Tj\^+s 


tr 




-(9+P+2) 


bN, 


i(uAr,«)dnj 


~ E f if *^0»^p('^i’'^2)^(a:i,a;2)uAr,«(dTi)j UAr,„(dr2)duj 

^ f e!^,l(AV’p)^z^iV,«(dT)duj 

^ 1 r=l In / 

2 771 / ^Tn+s r / \ \ 

EEEe / V’p(n, 'r 2 )^(aji, 2 ; 2 )„(' 1 - 2 , lAp)(dri)(dr 2 ) dn , 

p S 1 r=l ’ / 


Atv,! + AAr^2 + A7 v,3- 


We only treat the first term above and leave the remaining two to the reader. 


An,i = — 
ei 



1 

\^n\ 


E 

ieAjv 



E b)§(r)'Hp{Tj,U, 


Ti^u)'^{Xj,Xi) 


I m I ^Tjv+s 1 

" £1 5IAtvI 


2 

-{q+P+2),Ko,i,o 



2 

du 


where VAr^j,,- is the linear form introduced in (1,1.25^ for the choice of (ti^^,... ,tn,u)- By 
(j,1.29p . one obtains (writing w{Tk^u,Xi,u) instead of w{Tk,u,Ti,u, k^i, i-i)), 


Atv,! ^ 


Cm 

£i IAatI 



E w{Tk,u,Ti^u)'^{Xk.,Xi) 

keAff 


2 

dtt 


Cs 


£i IAatI" 


Y, df{xk,Xi)'i'{xi,Xi)E\ 




sup w{Tk^i 

u 


i)w{Tl 


,ui Xi^u) I ^ ) 

/ 


using Proposition l2.6l and Remark l.1..1l This term is lower or equal than any £2 > 0 provided 
s is sufficiently small. The treatment of the remaining terms can be done in the same way, 
using the continuity of the linear form Wat dehned in (j.1.26p . □ 

Theorem 3.23. The sequence of the laws of {TLn)n ^ 1 is tight in C{[0,T],WH|^’^°p^2)')■ 

Proof of Theorem \3.23[ Proposition 13.12l implies, as in the proof of Proposition 13.221 that 
the first item of the tightness criterion of Section 13.41 is satished for the process TLn- 
We verify now the second part of the tightness criterion. Since it is verihed for the 
process A4^ ^ (see Proposition 1,3.22p it suffices to prove it for the three terms (recall 
P-51P ) ^f'HN,sds, /o‘FAr,sds and /o*GAr,sds. Concerning ^f'HN,sds, if {Hp)p^i 
is a complete orthonormal system in '^'^Hp+ 2 '' have successively, for all £1 > 0 and 
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stopping time Tjv ^ T, 

'Tn+s 


K'HN,udu 

Tn+s 


cTn 

- ^*'HN,udu 


Tn 

Tn+s 

Tn 

Tn+s 


^:nMudu 


(II r 

^Hlfr 

s I r-'jv+s ^ 

<-E / Y. i'^N,u,^ui’pY 

£1 I JTn p j, 1 


-{q+P+2),Ko,Lo 


^ £1 


^:nN. 


(q+P+2),Ko,Lo 
du 


) 


“ »-(q+P+2),Ko,Lo 


du 


ef \Jtn 


I 'Hn,i 


-q,Ki,Ll 


dnj 


By Proposition 13.121 this expectation is finite. The same kind of estimates can be proven 
concerning the remaining terms using Proposition 13.16] and Remark l3.18[ The process 'Hat 
satisfies the tightness criterion, and hence, Theorem 13.231 is proven. □ 


3.6. Tightness of the fluctuation process r]]\f. Note that it is possible to retrieve the 
fluctuation process rjN from the two-particle process "Hat. Indeed, we easily see from (|2.9p 
that, for all function r /(d), 

{r]N,t , /) = , F) , (3.55) 

for the choice of 

F{0,uj,x,0,Cj,x) := d{x, xY f (9, Cj, x) (3.56) 

Unfortunately, deriving the tightness of r/Ar from the tightness of "Rw is not easy: the main 
difficulty is that, whatever the regularity of / may be, the function F in (|3.56p will always 
be at most a-Holder in the space variables (x,x), which is not a sufficient regularity for 
the test functions and the Sobolev embeddings used in the previous section. 

The strategy we follow is simply to adapt the proof of the tightness of "Rat to the set-up 
of r/Ar- Note that the proof becomes considerably easier since the definition of the process 
r/Ar (contrary to T-Ln) does not incorporate any spatial singularities: in the semimartingale 
decomposition (I2.11h of t/at, all the singularities in x are inside the term involving Hn- 
Due to the similarity between the proof of Theorem 13.241 and the arguments of Section 13.51 
we will only sketch the main lines of proof in this paragraph and leave the details to the 
reader. Recall the dehnitions of the weighted Sobolev spaces V in Section (3^ 


Theorem 3.24. For any q ^ P + 2, the sequence of the laws of (r/Ar)A>i is tight in 

C{[0,nV^XrPr2))- 

Proof of Theorem fS.Sfl Following the same procedure as in Proposition 13.121 it is easy to 
prove that sup^ ^ (|| m,t \\-q,Kuii) ^ H' W-gwiM stands for the norm given by 

(|3.7|) . A uniform bound about r]]\f similar to (|3.47l) . i.e. 


sup E sup 

i^N s: T 


'nN,t 


II 2 

II-(ij+P+2),ko,i.o 


< +00, 


(3.57) 


can be proven along the same lines as in the proof of Proposition 13.191 using the semi¬ 
martingale decomposition (12.lip of t/at- The continuity of the underlying linear operator 
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L[us] from o to V; 


q+P+2 


K.l,Ll 

q 


II llg,Ki,ti ^ C" II / \\q+p+ 2 ,K 0 ,L 0 ■ (3.58) 

can be proven as in Proposition 13.141 It is also immediate to derive bounds similar to 
Propositions 13.171 and I3.22] about the martingale process defined in (|2.13p (recall in 
particular Remark I3.18I) . that is 


sup E 

t ^ T 


M 


(v) 

N,t 


' ) 


^ c 


I^atI 


(3.59) 


The only unusual term in ()2.11l) is Jq {T~{-n,s , ‘^[/]) ds (recall the definition of 4>[/] in 
(j2.12p l. Using in particular the existence of a constant C > 0 such that !$[/]! ^ C || / ||g+p +2 
for any test function / e this gives that \{'Hn,s , ^[/])| ^ C || 'Hn,s \\-(q+p+2),Ko,io H ^ I 

and we can use the boundedness of the process "Rat in proven in Proposition [3T9] 

to conclude. The rest of the proof of Theorem I3.24I follows the same lines as in the proof 
of Theorem I3.23I □ 


4. Convergence of the fluctuations processes 

The rest of the paper is devoted to the identification of the limits as ^ oo of the couple 
in C([0,r],V © W) for some Sobolev spaces V and W defined in Section [T3l 
Obviously, since each component of this couple is tight, the couple is also tight. 

4.1. Identification of the initial value and the martingale part. 

Identification of the initial value. Note that at t = 0, the marginal ^o( d0, dw) on X x Y of 
the McKean-Vlasov process pq (II.8p is equal to C( d0) 0 //(dcu) (see (15.6h below). 

Proposition 4.1. For any q P + 2, under the assumptions made in Section lS.31 the 
sequence (r/Ar, 0 ) ”^^, 0 )^ of the fluctuation processes at t = 0 converges in law, as N ^ 00, 
^-lqlp+ 2 ) ® ^-(g+p+ 2 ) process {r]o,'Ho) defined by 

• if a < {r]o,'Ho) is a Gaussian process on ®^Yg+p+ 2 ) covariance, 

'^((|0>(^0)^=C^.(/l>/2)+Cl.,w(/l>52) + O,,p(/2,5l) + Cp(5l,I72), (4.1) 

where 

Crj{fij 2 )-= ^ f^{{fi{-r,x)f 2 {-,-,x), f,o)-{fi{-,-,x), (,o){f 2 {-,-,x), Co)}dx 

■= ^ f^[{f{-,-,x)[g^,fo]{-,-,x), ^ 0 ) - {f{-,-,x), f.o) {[ 9 '^,^o]{-,-,x), ^o)}dx 
Ch{9i,92) ■= ^ ^|([c/i4',Po] {■,-,x) [92'^,Po] {■,-,x ), Co) 

-{[ 9 i^Yo]{-,-,x), (o){[92^Yo]{-,-,x) , Co)}dx 

(4.2) 

where f\, f 2 and f (respectively 91 , 92 and g) are test functions on X x Y x § 
(respectively (X X Y X ), (/(■,■, x), Co) = /xxY^)? o( 40, dw) (recall (|2.1I1 ) 
and where [g^, vq] { 6 , ui, x) = [g'^, i 2 q ]2 ( 0 , u, x) (recall (12.3|) ). 


g+P+2,K;o!^0 
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if a>\, {rjOjT-Lo) is a deterministic process given by 
Vo = 0, 

{'Ho , g) ■= x(a) f 


/(XxY)2xS 


g{6,uj,x,6,Cj,x)^o{d6, da;)^o(d0, da;)da:, 


(4.3) 


for any bounded bounded function g with bounded derivatives and where x(«) 
is the constant defined in Lemma 


Proof of Proposition 14.11 is given in Section 15.5 


Identification of the martingale part. We identify in this paragraph the limit of the mar¬ 
tingale part M-n ■= ill the semi-martingale decomposition of the process 

{vn,Hn) (12.lip and (j2.14p . Here AIat is seen as a square integrable martingale in 
^-lq+P+ 2 )^ for any g ^ P + 2. 


■^^05^0 

-(q+P+2) 


-{q+P+2) 


Definition 4.2. Define the process M. = e 

lows: 

(1) If a < ^, M. is the Gaussian process with covariance, 

((), (^^)) E (/2)+, 

= (/i j2)+ 

where V^+p +2 ® ^TrPr 2 ’ e [0,T] and 


foi- 


(4.4) 


X X 1 f'S/\t n ^ 

Hif{fl,f2)-=-j^ j ^^df)^r)fl{T)dQ(r)f2{T)Uu{dT)du, 

/ \ ”1 f'Sfyt f' ID' 

(/>9) :=y E90M/(f)[((9g»5r)4',i/„]2(f)r'«(df)dR (4.5) 

1 r ^ 

lC^\gi,g2)-=- / XI [(^0»ffi)^>i"«]2('^)[(<9e»52)^,i"«]2(t^)i^«(df)du. 

Z JO J 


(2) If a> M. = () is equally 0. 

Proposition 4.3. For any qfi P + 2, the process M.n '■= (Ad^\ Ad^^) at ^ i is a martin¬ 
gale in ^T(^qlp^ 2 ) ®^-(g+p+ 2 ) Doob-Mcycr proccss iM.N)) with values in LO^ql’p+ 2 ® 

g^ven by 

+ iM'-FUgOiaV, 


In Eg, (A), and 

iMfuh)(h)=(Mfh , Mi;’*),, 

of rt m 

= 77^ E / T,dg^r)fi{0i^s,<Wi,Xi)dg(r)f2{0i,s,u;i,Xi)ds 
IAatI ieAjsr ^=1 


(4.7) 
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is the Doob-Meyer in associated to Doob- 

Meyer process associated to already calculated in (13.521) and 

(j2 m 

= rr^ E L T.^e(.-)fi^j,s)de(-)9{Ti,s,Tj,s)^{xi,Xj)ds 

IAatI ijeAjv-^° r=l 




i^A^r 


rt m 

ieAjv r=l 


where is given in (I3.54|) . Moreover, under the assumptions made in Section the 
process {Mn)n^ i converges in law in ® ^-(g+p+ 2 )’ o,s N ^ oo, to the process 

M given in Definition\4-S\ 


Proof of Proposition \4-3[ Estimates (BSI), dSD and (|4.8p are a straightforward conse¬ 
quence of the definitions of the martingale terms in (12.131) and in ()2.20l) 

and of the independence of the Brownian motions ,B^^^), i € Apf. It re¬ 

mains to prove the convergence of AIat to A4. First remark that the case a > | is a trivial 
consequence of the estimates (|3.49p . (|3.59p and the definition of oat in (ll.lOp . 

In the case of a < d, qn = \/N and (Al 7 v)Af^i is a sequence of uniformly square- 
integrable continuous martingales which is tight in V^ 2 (p+ 2 ) ® ^- 2 fp+ 2 ) • 
accumulation point. By arguments analogous to the ones continuously used in Section [3] 
(see in particular Proposition I3.12|) . it is easy to show that, first, defined in (13.541) 
vanishes to 0 as ^ oo and, second, that for all regular test function ( ^ ), one can make 
^ oo in the expression of ( g ) in (|4.6p and obtain for all t e [0,T] 

_ 1 p m 

((^))*(D=2 7o J E 

J f-t ^ fn 

^ / 'Z^e(r)f{r)[{dQ^r)g)'^,i^s]2{^)’^sidf)ds 

u X r=l 

1 r T^. n 

+ 2 70 7 E 


So, A4 is a square integrable martingale whose Doob-Meyer process is deterministic, given 
as above. So A4 is characterized by the Gaussian process with covariance given in ()4.4p . 
The convergence follows. □ 

Proposition 4.4. In the case of a < the processes (po^'Ho) introduced in Proposition \4. 1\ 
and Ai = ) defined in Definition \4-^ are independent. 

Proof of Proposition \4.4\ It is a straightforward consequence of the mutual independence 
oi {6ifi,Ui,Bi), i^ Km- □ 


4.2. Identification of the limit. So far, all the analysis has been made for general 
Sobolev indices {q,no,iQ). We now specify the adequate parameters: recall the definition 
of (n, t) in (12.261) and (|2.27p and define 

q:= P + 2 and (^ 0 ,^ 0 ) - (n + TT + t); 


(4.9) 
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SO that by definition of (ki, ti), ki = R and ii = l (recall (I2.26p and (j3.20p ). The conclusion 
of Section [3] for this choice of parameters is that the process {r]N,'HN)N is tight in 
C([0,T], ® ^- 2 (P+ 2 )) accumulation point {r],'H) belongs to 

the same space. In order to have a closed formula for (rj,T-L), we need to increase the 
regularity of the test functions, and hence, enlarge the space of distributions: we will give 
a characterization of in the larger space C([0,T], ®^l 3 (p+ 2 ))’ 


Proposition 4.5. Under the assumptions of Section \2.^ the process {r]N,'HN)N ^ i has 
convergent subsequences in ®^- 2 (P+ 2 )) accumulation point 

{r],'H) is a solution in ® ^- 3 (P+ 2 ) following system of coupled equations 


'ht = 'no + J^ L[i^sY'nsds + ^*nsds + m['^\ 
nt = no + ^fUs ds + , 


t e [0,T], 


(4.10) 


where (lyoj^o) is defined in Proposition \f. 1\ is given in Definition \4-S\ and 

L[i^s]* (respectively and ) is the dual of L[iys] defined in (|2.6p (respectively of ^ 
defined in ()2.12p and defined in p2.15p ). 


Proof of Proposition \4T^ Consider a subsequence (that we rename by {r]N,T~LN) for sim¬ 
plicity) converging in C([0,T], ®^- 2 ^p+ 2 )) We easily deduce from 

p3.47p and p3.57l) that 

e(sUP (|| Ht ||^2(P+2),«+7,i+7 + II ll-2(P+2),^+7,i+7)) < 

An application of (I3.34|) . (I3.37p and (I3.58|) in the case of g = 2(P + 2), kq = A and lq = l 
leads to 

II L[iys]f ll2(P+2),K+7,i+7 ^ C* II / ||3(p+2),K,t > 

• II ‘^[/] ll2(P+2),K+7,t+7 ^ ^ II / ll3(P+2),K,t ; (4-12) 

.11ll2(P+2),K+7,t+7 ^ ^ II 5 ll3(p+2),K,t • 

Estimate (14.lip together with (I4.12|] shows that /g* ds and Jq ^*'Hsds (respec¬ 
tively Jq P£f PLs ds) make sense as a Bochner integral in (I'ospectively in ^Y^^p^. 2 ))' 

Furthermore, Proposition 14.11 and 14.31 concerning the convergence of the initial value 
{ilNfl,T~l.N,o) and the martingale part (A4^\A4^^) and Proposition 13.161 are also valid 
for the choice of the parameters q = 2{P + 2), kq = k and lq = l. Moreover, for every 

ia) ^ ^3(P+2) ® ^3(P+2) ’ 

{VN,t, f) = {VN,0, f) + f {VN,S, L[’^s]f) ds+ f {nN,s,Hf]) ds + M^Y\f^ 

\ Pt pt ’ , 

, g) = {T-lNfi, g) + j {'Hn,s, ^sg) ds + j Fjsf,sgds + j GN,sgds + M^^lg. 

]4.13) 

For fixed ( ^ ) in ^ 3 (p+ 2 )®^ 3 (P+ 2 ) ’ define the continuous functional from C([0, T], V 7 g^p^ 2 ) ® 
W-f(p, 2 )) to m 2 by 

n C/ 7 d ■= I ~ -^0 ’ L[u,]f) ds - /o* {vs , $[/]) ds\ 

’ V {vt,g)- {vo , g) - Jq {vs , F^’sg) ds / ' 
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Since (0 in 


3(P+2) 


I w-’- 

'^'^3(P+2) 


^K+7,t+7 

2(P+2) 

,K+7,t+7 


i W; 


2{P+2) since is an accumula¬ 


tion point of {r]N,T-LN) in C([0,T], ®^- 2 (P+ 2 ))’ ^f,gi'nN,'HN) ^ 

as ->■ oo, which ends the proof of Proposition 14.51 □ 


4.3. Some uniqueness results. We derive in this section two uniqueness results that 
will be useful in the following. The first one concerns functional equations driven by the 
linear operator defined in (|2.16p . The second one addresses the same problem about 
the linear operator L[us] given by (12.6p . We treat in details the case of in Section IT^ 
The treatment of the operator T[Ps] being strictly identical, we only state the result in 
Section Ol and leave the proof to the reader. 


Uniqueness in linear functional equations driven by . Recall the definition of the 
operator in (I2.16p . For all 0 ^ s ^ ^ T and (r,f) e (X x Y x S)^, dehne as 

Xs,t{T,f) ■= (0s^t(r),0s,t(f)) the unique solution of the following stochastic differential 
equation in 


0 s , i ( t ) = ^ + v{t,&s^riT),uj)dr + Bt- Bs, 
0s,i(p) = ^ + v{t,&s^r{T),uj)dr + Bt- Bs, 


(4.14) 


for two independent Brownian motions B and B in (r, f) = {9,uj,x,9,oj,x) and where 


v{t, 0, w) := c(0, w) + [FT, ut] {Q,u). (4-15) 

We recall here that the quantity [FTjUi] is indeed independent of x (Remark [33]) . Define 
for any test function g 


U{t,s)g{T,f) = Epp(5r(Xst(r,f))). (4.16) 

Proposition 4.6. Under the assumptions of Section \2.^ for any functional R mC([0,r], W 
there is at most one solution in C([0,r], lo the equation 

St= £sds +f\sds, te[0,r], (^-l^) 

Moreover, one has the representation 

£t = £u{t,syRsds, (4.18) 

where U is given in (I4.16F 


Proof of Proposition 14.61 is given in Section 15.61 


Uniqueness for linear functional equations driven by L[ns]. We state here a result similar 
to Proposition 14.61 for the operator L[ps] defined in (12.6p . Since the proof is identical to 
Proposition 14.61 we only state the result. 

For all 0 ^ s ^ t ^ T and r = (6(,a;,x) e X x Y x S, define as Ysg{T) the unique solution 
of the following stochastic differential equation in X 

Ys,t{T) = 9 + u(t,i;,r(T),a;)dr + Rt-Rs, (4.19) 

for B Brownian motion in X and v defined as in ()4.15p and introduce the flow (for any 
test function /) 


P(t,s)/(r) = Ep(/(ni(r))). 


( 4 . 20 ) 
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Lemma 4.7. Under the assumptions of Section \2.2[ the operator Llvt] is continuous from 

'-^3P+8 '-^3(P+2) 


L[vt\f 


c'. 


0,2i. 

3{P+2) 


^ C II / 11^0^ _ , t e [0,t] 


\L[us]f-L[iyt]f\\ ^ 0 , 2 . ^CII/lUo,. |t- s|, e [0,r]. 

'-'3(P+2) '-"3P+8 

0,0 


For any j ^ 3P + 8, the operator V{t,s) is a linear operator from Cj’'^ to Cj''' such that 
II y(t, s)/11^0,1 ^ C II / 11^0,0 , 0 ^ s ^ ^ T, (4-21) 

^3 '^3 

I V{t, s)f- V{t, s')f I „o,i ^ C II / |Lo,o \/s' - s, O^s^s'^t^T. (4.22) 

'^j+i 

Moreover, the following backward Kolmogorov equation holds: for every f e 

4^(L'S)/-/ = ^ L[ut]V{t,r)fdr, in (4.23) 


One can deduce from Lemma 14.71 the proposition 

Proposition 4.8. Under the assumptions of Section \2.fA for any functional R inC{\d,T^,'Vz^p^^-^), 
there is at most one solution in C([0,T], equation 

St = £ L[ntYSsds +f\sds, te[ 0 ,T], ( 4 - 24 ) 

Moreover, one has the representation 

£t = £v{t,syRsds, mC°2°p^2)’ (4-25) 

where V is given in (I4.20p . 

4.4. Convergence when a<^. 


Representation of R in terms of r]. The purpose of this paragraph is to prove that the 
representation result (j2.inp of the two-particle process R in terms of the fluctuation process 
rj is indeed true when a < ^ (Section 12. ip . More precisely, we prove 

Proposition 4.9. Suppose a < 4. Under the assumptions of Section \2.2[ the process 

£N,tg ■■= {RN,t , g) - {gN,t Aut,^g)), g^ , t e [o, T] (4.26) 

converges in law to 0, in C([0,T], )’ as N ^ cx>. 

Remark 4.10. Let us admit for a moment Proposition \4-9[ This result (applied to the 
function g = 4>[/] for any r /(r) where is defined in (I2.12p ) combined with the 
analysis made in Section \4.2\ shows that each accumulation point r] of is actually a 
solution in uncoupled equation (12.361) . 

Proof of Proposition If-Hj The strategy is the following: first, write the semimartingale 
decomposition the process £n (based on the corresponding decompositions of pN and 
Rn), second, prove that any accumulation point is a solution to a linear PDE and third, 
prove that the unique solution to this equation is the trivial solution 0. 

First observe that the tightness of the processes Rn and rjN implies the tightness of 
the process £n itself in C([0,T],^ function (r,f) !->■ g(T,f). An 
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application of the semimartingale decomposition of (|2.11l) to a regular time-dependent 
test function (t,r) ft{T) gives 

{m,t, ft) = {VN,0 ’ /o) + ^ {m,s, L[us]fs) ds + ^ {nN,s , ^[/s]) ds + 

+ {'nN,s , dsfs) ds. (4.27) 

In the particular case of 

ft{f) = , (4.28) 


we have, using (HZl), 


, 4'(x,x) {^I\eg{T,f) + Veg{T,f) ■ \c{e,u:) + vt (r) 


and 


L[T^t\ft{f) = |i^i(dr), 4 '(x,x) QAg5r(T,f) + V0fi((r,f)- {c(0,£i}) + , vt ('^)}))- 

By definition of 4>[-] in (I2.12|] . we have also that 

^[ft\{T,f) = Veft{0,uj,x) ■r{6,u,e,u), 

= {I't, ^{■,x)VQg{-,T)) ■r{e,ui,e,Lj) = j^‘l-‘^'’g{T,f), 


where the definition of .jSf^ is given in (I2.17h . Consequently, one can rewrite (I4.27p as 
{VN,t, (ut, %)) = {VN,o, (Po , %)) + {r]N,s , 

(Wa-,. ,^Pg)ds + M^;ilft (4.29) 

Subtracting (|4.29l) to (I2.14p gives, for te [0,r], 

£N,tg = £N,og + i: £NA-^s^"^9)ds + i: FN,sg ds + i: GN,sgds + M^;^^g-M^;ilft. (4.30) 

We already know from Proposition 13.161 that F^^gg ds + G^^sg ds converges to 0 as 
A -> oo. We prove that it is also the case for the initial condition and the martingale 
part in (I4.30p . when a < |. First consider Sn^g = {'F-Nfi, g) - {VNfl, (^O; ^g))- Using 
ProDosition l4.ll we know that SnaO converges as A -> oo to {T-Lq , g)-{go , (i^o ; '^g))- Using 
the form of the covariance of (i/O; ^o) hr (14.2h . it is immediate to see that limTv^oo ^^7V,o = 0. 

We prove now a similar result for the martingale part in (|4.30p : we deduce from Propo¬ 
sition |43]that the martingale part in (I4.30p converges as N ^ oo to g-M^^'^ ft, where 
(_/\/((^), ) is given in Definition 14.21 It is then easy from the form of the covariance 

K-s,t in (14.5p and the definition of ft in (I4.28P to see that this process is equally 0. 

Following now the same procedure as in Section 14.21 we deduce from ()4.30l) that any 
accumulation point £ of £n in necessarily solves 

£t = f^\^g^^^)*£sds, te[0,r]. 

Then Proposition 14.91 follows directly from Proposition 14.61 □ 
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Convergence of {rj^) n^i- Theorem 12.71 is an easy consequence of the following uniqueness 
result. 


Proposition 4.11 (Uniqueness in law of the limit). Under the assumptions of Section \2.^ 
there is uniqueness in law (as well as pathwise uniqueness) in the limiting SPDE (j2.36|] . 


Proof of Proposition \4Ul\ We follow here the strategy used by Mitoma (see [32], Step 3, 
p. 352) who proved a similar uniqueness result in a slightly different context of weighted 
tempered distributions. Although the functional context is different, the strategy is simi¬ 
lar. Let t ^ r]{t) € C([0,r], ^ solution in dliM])- For simplicity, 

we write A4 in (|2.36p instead of in the following. Setting 


h{t) := Clr]sds = r]t-r]o-Mt (4.31) 

and differentiating this quantity with respect to t, one obtains that (almost surely w.r.t. 
the randomness) 

-^h{t) = C)h{t)+r]o +Mt (4.32) 

at 

Define for all 0 ^ t ^ T, all test function / and {9,uj,x) 


Ar[t't]/(6',a;,x) := {ut, V0/(-)' r(-,-,6',w)T(-,x)), 

so that the linear operator Ct in ()2.37p 

Ct = L[ut]+K[ut]. 


(4.33) 


(4.34) 


Focus for a moment on the regularity of the linear operator K in (I4.33p : for all / ^ 1, 


there exists a constant C > 0 such that for all / e for all u e (0,r] 

1 




■'i' ) ’ 


(4.35) 


where ao e [0, |) is dehned in Proposition [221 Indeed, using the decomposition vti d9, dn;, dx) 
pt{9, Lo) d9fr{ did) dx (recall Proposition 12.5|] . one obtains by integration by parts 

A'[^'«]/(6',w,x) = - f divQ{pt(9, Q)T(9, uj, 9, Ld)) /(0,a;,x)T(x,x)di9/i(da;)dx, 

JXxYxS 

= - [ divg (pt{9,Ld)T(9, Cd, 9, id)) f (9,0, X - z)p{z) d9iJ,(dCd) dz, 
JXxYxS 

where p is defined in (I3.36p . Hence, for all / ^ 1, 


II K[iyu\f Lo.o = 

Y, sup I r dwg{pt{9,id)DQ^D^^T{9,id,9,id))D^^f{9,id,x-z)p{z)d9p{did)dz 
|fc| ^ie,u},x IXXxYxS 

11/11^0,. Y f \dive{Pt{^^d))Dg^D^^T{9,u>,9,uj))\{l + \id\‘')d9p{di:d), 

‘ |A:| sj Z -'XxY 
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where the sum above is taken over all multi-indices k = {ki,k 2 ,k 3 ) such that |A:| = |A:i| + 
1 ^ 2 ! + 1 ^ 3 ! ^ b Consequently, using (I2.3.sp and (I2.34j> . 

II ^[^u]f lUo.o ^ 

^C||/|LoJ —+ ^x|, 

\m“o ii“o+^ / 

by the assumptions made on fj. and T (recall Section [221) • Using the semigroup V{s,t) 
associated to L[ui], one has 

-^V* {t,u)h{u) = -V* {t,u)L[iyuY h{u) + V* {t,u){C*^h{u) +% + 
du 

= V*{t,u){K*[vu\h{u) +m + Mu)- 


Hence, h is also solution to 

h{t) = f V*{t,u){K*[i 2 u]h{u)+r]o +Mu)du. 
Jo 


(4.36) 


The main point of the proof is to see that h solution of (|4.36p can be approximated by 
the converging sequence {hn)n^i defined recursively as follows 


{hi{t) = jQV*{t,u){r]o + Mu)du, 

|/in(t) =fQV*{t,u){K*[uu]hn-i{u)+r]o + Mu)du, n ^ 2. 


(4.37) 


Indeed, by the boundedness of the semigroup V{t,u) (14.211) and by the boundedness of 
K[i/u\ ()4.35p . we obtain that for all 0 < tt < t < T, for all / e , for all h e C^^p^ 2 ) 

\hK[i/u]V{t,u)f\ II^IUo’O II-^[t'n]U(f,u)/|Lo,o 


-2(P+2) 


^ (7 II 11^0,0 

^-2(P+2) 


^ (7II |Lo,o 

‘^-2(P4-2) 


\m" 0 li«0+2/ 

\tt“0 y^»0+2 ) 


"2(P+2) 


y{t,u)f Lo,. 


■^2(P-t-2) 


/ lUO’O 


"^2(P+2) 


Thus, the sequence {hn)n ^ i defined in (14.371) satisfies, for all n ^ 2 


^'n+l(0 ^n(0 




-2(P+2) 


rf— 

Jo tt“ 0+2 / 


^n(ll) l(ll) 11^0,0 


du. 


''-2(P+2) 


Choose now r > 1 such that 1 < r < —Ur and r* such that - + U = i. Then, by Holder in- 

oo+l r r* ’ 

equality, one obtains that || hn+i{t) - hn{t) ||U,o ^ C C || hn{u) - hn-i{u) ||()o,o dtt. 

^-2(P+2) -2(P+2) 

By an immediate recursion, for all A; ^ 1, || hn+i+kY) - hn+kY) ll^o.o ^ so that 


{hn)n^ 1 is a Cauchy sequence in C([0,r],C°2°p^2)) converges to h, solution of 

()4.36p . Turning back to rj (recall ()4.3ip ) and writing R{t) := rjQ + Mt, we obtain that r] is 


-2(P+2) 
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uniquely written as 

r]t= lim \ r]o + Mt+ f V* R{ti) dti 

n^oo Jo 

V^tRi)K[ut,rV^t,t2)jC-t,R{t2)dt2dh + ... 

X t rt\ rtn-1 

■■■ V*{t,ti)K[ut^Y--V*{t,tn)^„R{tn)dtn... dt2dtiy. 

(4.38) 

This proves pathwise uniqueness. But if one chooses another solution fj defined on another 
probability space, with initial condition fjo and noise A4 with the same law as {r]o,M), we 
obtain the same expression as above with R replaced by R{t) = fjo + Mt- Since R and R 
have then the same law, uniqueness in law in (12.361) follows from (|4.38l) . Proposition 14.111 
is proven. □ 


4.5. Convergence when a > |. Theorem 12.81 will be proven once we have established 
uniqueness of a solution to (|2.38p . The main issue here is that the linear operator 
in (I2.15P is not the generator of a diffusion, due to the nonstandard integro-differential 
term in (12.151) . In particular, the existence of a semi-group solving a backward 

Kolmogorov equation similar to ()5.52p is unclear in the case of ^g. Thus, we restrict here 
to the case where the state space X is compact and c bounded (see Section 12.41) . 


and\2.4[ there exists at most 


Proposition 4.12. Under the assumptions of Sections] 
one solution (r/,^) in ®'^- 2 (P+ 2 )) (|2.38p . 

Proof of Proposition \4-t2[ The supplementary assumptions made in the beginning at the 
paragraph ensure that the weighted norms introduced in Section 13.31 are now equivalent 
to the usual Sobolev norms without weights. The purpose here is to use the regularizing 
properties of the heat kernel in those spaces. 

Let be two solutions of (j2.38p in C([0,T],® 

is in C([0,T], solves 


^- 2 CP+ 2 ))‘ Then the difference PL 


in W-’- 

m vv_3^p^2) 

PLt = fJi^sYnsds. (4.39) 

In particular, it is a solution to the following equation (with unknown 6 e C([0, T], )) 




where 


Rg:= (^J^'g-^AgJi)|*ng, se[0,r], 
Denote by (0,9) g^{9,9) the heat kernel on X^ and by 


(4.40) 

(4.41) 

(4.42) 


G{s,t)ip{T,f) := gt-g * (p{T,f) 
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the corresponding semi-group (note that the convolution only concerns the variables 
{9,9)). It is standard to prove a backward Kolmogorov equation similar to (15.521) con¬ 
cerning G: 

G{t,s)g-g = £ g{G{t,r)g) dr, in C^p+ 2 )- (4.43) 

By the same procedure as in Proposition 14.61 one obtains that (I4.4U|) has a unique solution 
(which is necessarily Ti) that satisfies, 

nt = f^'G{t,s)*Rsds, inC_3(p+2), (4.44) 

that is, for all g in C 3 (p+ 2 )) 

(^ 1 , 5 ) = ^ {^'Hs,{^.Sfs-^AggjG{t,s)g'^ ds. (4.45) 

Note also that the coefficients in do not depend on the space variables {x,x) (by 
rotational invariance), and that 4>[/] in (12.121) only depends on x, not on x. Hence, with 
no loss of generality, we can consider test functions g depending only on {9,uj,x,9,Cj). 
Using the regularizing properties of the heat kernel, we deduce from (j4.44p that, for some 
constant C > 0 and some /3 e (0, 4), 

II ll-3(P+2) ^ X ll-3(P+2) 

Let us admit for a moment ()4.46p . Then, following the same procedure as in [13], Step 3, 
p. 763, we obtain directly the uniqueness for T-L. Consequently, fj := gG') _ solves 

fit= [ LlusYfisds, te[0,T]. (4.47) 

Jo 

Applying Proposition 14.81 we obtain that fj = 0 and Proposition 14.121 follows. Hence, it 
remains to prove (I4.46p . One has for any differential operators Di and D 2 of order qi and 
q 2 such that 3(P + 2), 

\\Di{S/ 0 G{t,s)g) ■ {D 2 c)\\ 1 = f \Di{VeG{t,s)g) ■ {D 2 c)f drdf, 

J(XxYxZ)2 

< drdf < ^ I D,g ||^ . 


The term Veg{x,f) ■ [FT, (p) in (I2.15h can be treated in the same way, using the 
boundedness of P and its derivatives. For the last term in (j2.15p . for the same reasons, it 
suffices to estimate, for any differential operator D of order smaller than 3(P + 2), 

f 4'(-,x)V0ff(-,r))f d6'da;da: = J Vgg{9,u!, z, 9,uj)ps{9,uj) d9g{ dw) dz 

X - z ■, •, z, 0,0;), dz| dxd0do;, 

" f i f \P* {x)dx] d9doo, 

JXxY \~/S / 


d0 do; dx, 
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where we recall the definition of p in (I3.36|) . By Young’s inequality, 
f Hi's, d6»da;dx ^ \\p\\l 


= \\P 


€ 


|?X 

|?X 

-/ 

3 Ja 



Vgg{6,uj, X, 6,uj)ps{9,uj) d9p{ dw) 

2 


dx I d9 dw, 


^ ll/o|li / \V§g{9,u},x,9,uj)\ ps{9,u!)d9p{dio)d9diodx, 

(XxY)2xS 

2 

^ ^ ^ 05 ( 0 ,a),x,6*,a;)|^ (1 + |a)|‘')d0d6*da;/r(dw)dx, 


s“o J(XxY)2x; 

where we used the a priori estimate (12.3311 on the density p. As in (14.441) . injecting 
DrG{t, s)g = G{t, s)DTg into the previous estimate, we obtain 

J \{vs , ^{■,x)VgG{t,s)Drg{-,T))H d9dujdx 

^ X X, 0,a;)|^ (1 + |w|'') d9 d9 dLjjp{ do;) dx. 


II Pi 

s“o 


Another application of Young’s inequality and usual estimate on the heat kernel lead to 

G 


J \{us, ^{■,x)VgG{t,s)Drg{-,T))f dPdwdx 

The inequality (14.4611 follows by density. 




s"o(t - s) 


n l|2 

^rg II2 • 


□ 


5. Proofs 


5.1. Proof of Propositions [272] and r2.31 We prove in this paragraph the semimartingale 
decompositions of processes pN and T-Ln of Section ITTI 

Proof of Proposition \2.2[ Writing Ito’s formula for (HI) gives 

't I 1 


{I'Np , f) = {l'N,0 , f) + , -Ae/j 

+ X {i^N,s,^ef ■ [c + [T^,m,s ]) ds +-^M^j^lf (5.1) 

Combining the McKean-Vlasov equation (ll.7jl and (15.ip , one obtains 
{VN,t, f) = {PN,0 ’ /) + X 

+ aivX {{i'N,s,yef-[^'^,’^N,s )Ve /-ds + , 

= iVNfl > /) + X ■ c| ds + M^^lf 

+ X [aNi^l'N,s-I's,yef ■[^^,1's '^ + aN{^l'N,s,y9f-[^^,l'N,s-l^s]))ds, 

= iVNfi ’ /) + X }) 

+ aN{^’^N,s,yef ■[^^,i'N,s-i^s |ds+A4^|/. 


'=y^N,t 
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The remaining term ujv can be computed as: 

/ 1 


UN,t = f 

Jo 


Oat 


2 S V0/(Tj^s) ■ Wj, Xj) 

\-^ n \ ijeAjv 


E f yefin,s)-m,. 

I^A^I ieA„ 


,Ui,6,oj)^ixi,x)i^sidf) ds 


= £ {nN,s,Hf])ds. 
Proposition 12.21 follows. 


□ 


Proof of Proposition \2.J[ Applying Ito’s Formula to (|2.4I) (for any regular two-variable 
function (r, f) i-s- (^(r, f) = gi 6 ,uj,x, 6 ,u},x)), one obtains, for all t e [0,T] 

1 r* 

9in,tiTj,t) — 9in,OiTj,o') 2 J g Tjjs) ds 

+ yegin,s,Tj,s) ■ {ci 6 i^s,u}i) + [rT,i/Ar,s](ri,s)} ds 


^§g{n,s,rj,s) ■{c{9j^s,^^j) + [r^',UAr,s](Tj-^)} ds 


(5.2) 


• dB 


y§gin,s,Tj,s) 

Applying Ito’s Formula to the process 1f ^(xj, x)( 5 :(rj^t, f)t't( df) (recall (Il.Tp i: 

J'i>ixi,x)giTi,t,f)utidf) = J T(xi,x)5(ri,o,f)uo(df) 

'^\lo I 

^ fo I ■ |c(6»i,s,Wi) + [rT,UAr,s] (Ti,s)| ds (5.3) 

^ fo I i^i^ sgin,s, T)nsi df) ■ dBi^s 

^ fo I + [r^',Us] (T)}us(df)ds 

Combining (|5.2|) and (|5.3p . we have 

i'HN,t,g) = {'Hn,o, 9) + {^'Hn,s , ^Agg5(r,f) + V0C/(r,f) ■ 0(0,0;) + VgC/(r,f) -0(0,65)1 ds 

+ {'HN,s,^9giT,f)-[T'Il,UN,s]iT)) ds + M^j^^g 

f \ f ^ixi,x)VggiTi,s,f) ■[T^,i^n,s] (f )i^N,si df) 
,= A -'0 \ J 


aN 

I^nI 


ieA N 


J 'itixi,x)VQgin,s,r) • [FT,!/^] (f)us(df)| 

(5.4) 
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We treat the last term of (j5.4p apart: it can be written as the sum of Ajsi and Bjsi where 


dN 


icA., ~'0 '' 


^N,t-= ,2 ^ 

\^n\ ijeAjv 

= [ E E ^§9iB,s,rj,s)^{xi,Xj)] ■ Y ^{dj,s,u;j,Oi,s,uJi)^{xj,xi) 

9o IAatI i-IAArl J 

-fnh. ,u]j,9,Cj)'^{xj,x)i's{ df)| ds, 

X t _ 

{'Hn,s, {i'N,s, Vggi-,T)'if{-,x)) ■r{e,uj,e,cj)) ds, 

and Bj\f is given by 

yiAivIjeAjv 

- J 'i/{xi,x)VQg{Ti^s,T) ■ (f)i/^(df)j ds, 

"X ds. 

Rewriting (I5.4p in term of A]\f, B]\f, Fn and Gn gives Proposition 12.31 □ 

5.2. Proof of Proposition [27^ A solution to (II.7p is provided by the nonlinear measure 
(■^i)te[o,T] introduced in Proposition 12.41 Indeed, applying Ito’s formula to any function 
{9,oj,x) ^ f{9,uj,x) that is w.r.t. 9, one obtains: 

f{9^’'^,ui,x) = f{9o,u:,x) + ^J^ ds + ^ Ve/■ c(0f’‘^,a;) ds 

+X* ^ X* 

Taking the expectation in pS.Sp gives that 1Xt is a weak solution to pi.70 . The uniqueness 
in pi.7p is a consequence of [28], Proposition 2.19. 

We now turn to the proof of the decomposition (I2.3ip : assume for a moment that a 
decomposition like (I2.3ip exists. Necessarily ^ must satisfy 

^o( d0, dw) = C( d0)/x( do;) (5-6) 

and for all t ^ 0 and regular test function {9,uj) h{9,u) independent of x 

dt (6 , h) = + Veh ■ c| + {ut , Veh ■ [p^', ut |. (5.7) 

Since by assumption, vt{ d0, dw, dx) = ^t{ d0,a;) dx, 

[P'kjt'i {9,u},x) = [ r(0,u;, 0,d;)'I'(x, x)^t( d0,w) dx, 

L -I ^XxYxS 

= ^T{9,uj,9,oj)Ct{d9,oj)^ ^^^'(x,x)dxj . (5.8) 

As already noticed in the proof of Lemma [3~il the integral /g 'I'(x, x) dx = 2 
is independent of x. Consequently, m becomes 

dt {f,t, h) = ^Aeh + Veh ■ c| + , YT^^eh ■ [p, (,t |. 


(5.9) 
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Equation (|5.9|) endowed with ^o(d6*, dw) = C{d6)fj,(duj) is a standard McKean-Vlasov 
equation with regular coefficients, and hence, admits a unique solution (see for example 
m or [331? Lemma 10). Defining Xt(d9, do;, dx) := Ct(d9^ da;)dx and observing that 


2^ 

l-a 




r^',v 


one readily sees that 1Ai is a solution of (jl.7p . By uniqueness in 


(jl.7p . A = u and (12.311) follows. 

We focus now on the regularity of solution of (|5.9I) . We follow here the strategy 
developed in m, Proposition 7.1: fix T > 0, a; e Supp(//) and 1the unique solution in 
C([0,r], Ali(X X Y)) to (j5.9h . Define At{0,cj) := c(9,uj) + {9 ,lo) and consider 


the linear equation 

dtpt{0,uj) = ^Apt{0,uj)-diY0(j)t{9,u})At{9,u})Y 
such that for p-a.e. uj, for all / e C(X), 

f^f{9)po{9,u;)d9 = f^f{9)C{d9). 


(5.10) 


(5.11) 


For hxed uj e Supp(/r), A.(-,uj) is continuous in t and regular in 9, by assumption on P. 
Suppose for a moment that we have found a solution pt{9,uj) to (I5.10p - (I5.11I) such that for 
/x-a.e. u!, pt{-,uj) is strictly positive on (0,r] x X. In particular for such a solution p, the 
quantity f^pt{9,uj) d9 is conserved for t > 0, so that pt{-,uj) is indeed a probability density 
for all t > 0. Then both probability measures d9, dw) and pt{9,uj) d9p{ dcu) solve 


f f{9,uj)iut{d9,duj)= f f{9,uj)C{d9)fi{duj) + ^ f f A 0 f{ 9 ,uj)us{d 9 , duj) ds 

JXxY JXxY 2 Jo JXxY 

+ f f V0f(9,uj)As(9,uj)iys(d9,duj)ds. (5.12) 

Jo JXxY 

By m or [33], Lemma 10, uniqueness in (15.911 is precisely a consequence of uniqueness 
in (I5.12p . Hence, by uniqueness in (15.121) . ^t{d9, dw) = pt{9,uj) d9p{du}), which (12.321) . 
So it suffices to exhibit a weak solution pt{9,uj) to (15.101) such that (|5.11l) is satished. 
This fact can be deduced from standard results for uniform parabolic PDEs (see 13 El 
for precise definitions). In particular, a standard result (see [IT], Theorem 2, page 251 
or [3|, Section 7) states that (I5.10p admits a fundamental solution G{9,t]9', s,uj) {t > s), 
satisfying, for all j € {0,1}, for any differential operator Dj ^0 in 9 of order j 


\Dj^eG{9,t;9',s,uj)\ ^ 


__.expf-J^) 

{C2{uj){t - S))"^ V C2{uj){t-S)l' 


(5.13) 


for positive constants ci and C 2 that only depend on the bound and the modulus of 
continuity of At (see Theorem 1, page 241). In particular, by the assumptions made 
in Section [221 one can suppose without loss of generality that for all f = 1,2, Ci(uj) = 
Q;i(l V |cu|^), for some at > 0. Thanks to Theorem 3, page 256 in [17j . the following 
expression oi pt{9, uj) 

pt{9,uj) = J^G{9,f,9',0,uj)C{9')d9' (5.14) 

defines a classical solution of (I5.10p on (0, T] xS such that (15.lip is satisfied. The positivity 
and boundedness oipt{-,uj) for f > 0 is a consequence of [1], Theorem 7, page 661. Applying 
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(j5.13p with j = 0 gives 

0^pt{0,uj) ^ m f ew{--^- 7 ^\c{G')de', 

JX \ C2{0J)t f 


€ 


ic2iuj)t) 2 
ciito) 


, , . , , exp, 

{C2{0j)t) 2 \JX \ C2{u!)t 






{C2{i0)t)- 


{C2{oj)t) 2 9^ 


where in (|5.15p ^ ^ = 1 (recall ()2.281) 1. Similarly, using again (I5.13[l . one has 


|divepf(6»,a;)| ^ C- 


ci(a;) 


l,m/i 1\ 1\' 


(5.15) 


(5.16) 


The parameter p in (I5.16P is such that ^(1-^) < \ (recall p2.28p l. By definition of the q, 
i = 1,2 and by (12.291) . the estimates p2.33p and (12.341) follow. Proposition [23] is proven. □ 


5.3. Proof of Proposition [3T6l Since the initial condition and the Brownian motion are 
the same for 6i and 9i, one has: for all i e A^v, 

^ 32 ^^ {ei^s-di,s)-{c{Oi^s,u}i)-c{ei^s,^i))ds^ 

+ 32^J^ (6'i,s-^j,s) ■ -^T^,i2s (0j,s,Wj,Xj)) dsj . (5.17) 


Using the one-sided Lipschitz condition on c, one can bound the first term in (I5.17h by 
32L^ (/o' \9i^s - Oi/ ds)" ^ 32L^T^ // ds ^ 32L^T^ // max^.A^ sup„ ^ - 0^,/ ds. 

Now, introducing the empirical measure of the nonlinear processes (0i,..., 0Ar) 

= p/i E Sr,,, N»1 (5.18) 


jeAjv 


the last term in p5.17l) can estimated above by C{A]\f + Bn + Cn) {C being a numerical 
constant) where 


— \0i,s-di,s 

B „(/ \e, 

Cn '■= \ 0 i,s - I 


-' 2,5 


PT 

PT 

PT 


1^N,S {9i^s,i^i,Xi) - ^T^,Un,s {9i,s,^i,Xi) dsj , 

^N,s~^N,s (,0i,sj Xi') ds^ , 

VN,s-T^s {9i,s,^^i,Xi) dsj . 


(5.19) 
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Concerning A]\{ and B]\j, we have, by the Lipschitz continuity of F, for some constant 

C>0, 

An ^ [£‘h,-k.\ J\m.. ,uji,9,uj) dw, dx)dsj , 

-t 


€ II rr.. 


Lip 


< 


CT^ II r 11^. f max sup \9iu - 9i^uf ds, (by ([33])), 

Jo leAjv u<s 


and, 


Bn ~ ( lA I f |^*,s ^i,s| Wj, s, cUj) Wj, Wj)| ^(xj, Xj) ds I , 

\ -JO j 


^ II r 


lll,ip(|A^| XI '^iXi,Xj) \9i^s %,s| dsj , 


^ II r 




I Lip I 


max sup \9i^u - ds 

i^^N Ft < S 


JAnI 

CT^ II r 11^. f max sup - 9iJ^ ds, (recall ([33])). 
^ Jo ieAiv u ^ s 

It remains to study the term Cn (I5.19P . Since we obviously have 


)' 


_ a 7^3 r~t 

Cn^— K«-0*,srds + — 

2 Jo ' ' 2 Jo 

r*‘ I - i8 r^w - ^ - 

^ FT / maxsup ds+— r^',UAr s-i/s (9i s,uji,Xi) 

2 Jo ieAjvnCs 2 Jo IL J ’ 


T^,Un,s-’^s {9i^s,<Ji,Xi) 

r* 


ds 


we concentrate on the term cn ■= 
positi 

(In ■= 


r^r 




ds, (5.20) 
9i^s,Ui,Xi) in (|5.2UI) . Using the decom¬ 


position u{ d9, du, dx) = (,{d9, dn;)dx (recall (I2.3ip l. we have cn ^ Jn + &N where 
1 


(5.21) 


Cat := 


y ^{xi,Xj){T{9i^s,0Ji,9j^s,^Jj) - / r(6'i,s,a;i,6',a;)^s(d6', dcj)] , 
jeAjv \ J / 

^ H f ^{0i,s,^Ji,0,ij)^s{d9, doj)di{xi,Xj) - f T{9i^s,uii,9,ui)^s{d9, du})^{xi,x)dx 

e\.rJ J 

(5.22) 


I^Af| 




We treat first the term d n ([5.211) . For simplicity, we use the notation Ti j := r{9i,uji,9j, loj )- 
im ,uji,9,u})£s{d9, dcj). Note that, by independence of the 9i and by definition of Tij, 
for any J ^ fci, /c 2 , /cs and i J 

■ Ti^ki){Ti^k2 ■ = E(E((Tj^j ■ Ti^ki)iTi,k2 ' Ti^k3)\Gr,r * J)), 



\ 


E 

TiM ■^{Ti,j\9r,r t l) 



= 0. 


=0 
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Consequently, 


e(4) = e 


jz^i 


, E fn ^{\TiMf \TiMf {TiM^ \T. 

I^ivl ki,k2MM*i^l=^ ' 

1 


i,kf) 


^ (n E(|T,,fcj' \Ti^ksf 

:2,Ai3,fc4 \/=l / 


\^n\ ki,k2,ks,k4 
distincts 

^ |r,/) 

\^n\ j,k,l 

distincts 

^ ^(Xi, Xk)^^{xi, xi)^E{\Ti,kf \T^,l\^) + E E \Ti,kf ^^^ix^Xk) 


lAivr 


k^l 


^ (CII r I 


1 


E + 


JAatI fci,fc2,fc3,fc4«=l 

distincts 

IAatI k*i 


IAatI k 

-E— Y, ^{Xi,Xjf4!{Xi,Xkf^{Xi,XlY 

IAatI j,fc,z 

distincts 


1 


IAatI k 


^^'(Xj,XA:) 


8a 


(5.23) 


where we used that iTjjl ^ 2 || F and where C is a numerical constant. An application 
of Lemma [32] on ()5.23p leads to an upper bound of E(d^) of the form ^ if 0 ^ a < |, 
(respectively if | < a < 1), for some constant C. It remains to treat the term bat 

in (|5.22p : for simplicity, we use the notation ni^s ■= J ^{di,s,‘Xi,9,u})^s{(id, dcu), so that 


Bn = 




E f T^i,s{'^{xi,Xj)-^{xi,x)) dx ^ ||r||^ Y f \^{xi,Xj)-'d>{xi,x)\ dx, 

eAjv jeAjv 


jeAjv ^ 

g||rL 

Ni-a 




Taking the expectation in (I5.17P and applying Gronwall’s lemma ends the proof of Propo¬ 
sition 13.61 □ 


5.4. Proof of Proposition I3.12L Since we work in this proof with the same Sobolev 
space (and its dual), we write in this proof || • instead of || • ||g for simplicity 

of notations. Introducing the nonlinear process fi^t = {9i,t,^i,Xi) fProposition 13.61) . one 
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can decompose ■, d) into the sum of three terms: 


aN 


IAatI ijeAjv 


(5.24) 


'JAivI jjeAjv IAa^L>a..D 


|An| 


zeAjv 


K^nAs) ■=E f '^{Xi,x) {g{T^,f) - g{Ti^t,f))’^t{df), 

IAa^I jeAjv 


(5.25) 

(5.26) 


Intuitively, the boundedness of In and /Cat in this decomposition comes from Proposi¬ 
tion [321 each particle 6i is well approximated by its corresponding nonlinear process 
The most technical part of the proof concerns the evaluation of In and relies on cancella¬ 
tion arguments in the expression of In{9A, using the independence of the (0i)ieAjv The 
main point of the proof is that we need an estimate that is uniform in g. To do so, we use 
Hibertian techniques similar to miET]. 

Concerning In Ad)-: 


^nA9? 


‘'N 


|Aiv| i,j,k,l€A 


E '^{Xi,Xj)'i>{Xk,Xl) {g{Ti^t,Tj,t) - 9{Xi,t,TjA) {9{Tk,t,Tl,t) - 9{Tk,t,TlA) 




E ^{xuXj)'^{xk:Xi)TZi^jA9)T^k,iA9)^ 


|Aiv| ij,k,l€AN 

where we have written, for simplicity (’^scall ()3.2ip b Applying the 

previous equality to an orthonormal system {gp)p i in the Hilbert space and using 

Parseval’s identity, one obtains 


\^N,t\\\=\ E \'^nA9p)Y 

\p5= 1 


77^ E '^{xi,Xj)^{xk,xi) E T^i,jA9p)T^k,iA9p) 

[AatI i,j,k,UAN I 


< 




lATvr 




E T(xi,Xj)T(xfc,xO II ||7^fc,^,t|| , 




,, ,8 E 

[Aivl /=! 1=1 


Jlvj4 


By Proposition 13.101 and Cauchy-Schwarz inequality. 


E 


(n 


7^. 




< cnE(x*,,,,)’ nE((!«,,,-+ \e „-< A, 
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where Xi,j,t is a shortcut for (recall (j3.27p l and where we used ([37 

and Proposition 12.61 Consequently, using (|3.3|) . 


sup E IIX, 


1 ^ 7V,t ^ T 


'-N,t 


) ^ sup 
1 AT \ 


C 


E n E <+0°- 


I^A^I i \. .jideAjv ^-1 \\-^n\ 

j<^^N 


We now turn to the estimation of JM,t in (|5.25l) : for all ^ 1, t e [0,T] 

JN,t(5) = 7^ E E ^ixi,Xj)g{fi^t,fj,t) - f ^{xi,x)g{fi^t,f)i't{df)], 

|i^Af| ieAjv \ ieAjv / 

= 77^ E ^ixi,Xj)ig{fi,t,fj,t) - g{fi,t,0,u,Xj)^t{d6, du)] 

IAatI ijeAiv \ J ) 

+ E 1 E ^ixi,Xj) f g{fi,t,9,oj,Xj)(tid9, did) - f ^{xi,x)g{%t,f)ut 

I^A^I ieAjv \ ieAjv 

^=4'i(5) + 4'i(5)- (5.27) 

The linear form in (|5.27l) captures the random fluctuations of the nonlinear processes 

around their mean value (for fixed positions), whereas ^ captures the variations of the 
positions of the particles (for a fixed randomness). The first term in (I5.27P can be simply 

written as J^j-ig) = Eij.Ajv ^ixi,Xj)Aij^t{g), where 

AjAd) ■= g{n,t,fj,t) - J g{fi^t,9,ui,Xj)Ct{d9, dw) = g{fi^t,fj,t)-'E{g{fi,t,fj,t)\fr,t,r* j ). 

(5.28) 

Using again a complete orthonormal system {gp)p^ i in 
/ \ 2 


J 


( 1 ) 


N,t 


‘ I E 


l r ^ 1 


'‘N 


E E '^iXi,Xj)'dl{xk,Xi)Aij^t{gr)Ak,l,t{gr)] , 


VIAatI r 3 s li,j,fc,Z£Aiv 
q4 4 

-Pw E E UHxn, Xji )Aq ji(g'r)Ajj t{gr)Ai^ )Tj^ ,ji,t{gs )• 

IAatI r,s Js Z=1 

j'lvir 

(5.29) 


Taking the expectation in (j5.29l) . we need to estimate the quantity 

E{Il, I 2 , Ist I 4 ) ■- ^ {^ii,ji,t{gr')Ai2j2,t{9r')-^i3,j3,t{gs')-Ai^J^^t{gs')') 1 (5.30) 

for every sets of couples /; := (f^, j)), / = 1,..., 4. One crucial observation here is that for 
all j € Atv, 

^{^i,jA9)\xr,t,r * j) = 0, (5.31) 

by definition of the Aij in ()5.28l) . ()5.3ip leads to cancellations in (I5.30p . We give here 
a simple combinatoric argument in order to derive the relevant contributions to the es¬ 
timation of the expectation of (I5.29p : associate a vertex to each {ii,ji} (/ = 1,...,4) 
and draw an edge between two vertices labelled k and I {k,l = 1,...,4) if and only if 
Ahji} ^ {ik,jk} * 0 (see Figure [I]). Note that when a vertex (labelled for example by 
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Figure 1. The only nontrivial contributions to ()5.30l) corresponds to in¬ 
dices given by the above diagrams. 


is isolated (that is n {i 2 j* 35 * 45 j 2 ) jsj jd} = 0), by independence of the 

random variables fj, one obtains that 


F'(/l, I 21 IZl ^(^*2 J4 


E(A 


UJil'r-) 


r*ji)) = 0. 


Using this last observation, it is easy to enumerate by inspection all the possible nontrivial 
contributions to the expectation of (|5.29p : they are listed in Figured) We only give here 
a detailed analysis of the two extreme cases (i.e. the largest (A) and smallest (G)) and 
leave the remaining cases to the reader. The largest contribution of E{Ii, 12 , 13 , 14 ,) in 


El 


J 


( 1 ) 


N,t 


-Q/ 


1 corresponds to the case (A) in Figured) where the set of indices {Ii, I 2 , 13, 14 } 


is separated into {/i, 12 } and {I 3 , 14 } with the property that {ii, ji, ^ 2 )^ 2 } {^ 3 ) J 3 U 4 )^ 4 } = 
0. By independence of the random variables {fi)i^A^, the contribution of this term to 




is then 


E 77 ^ E '^ixi,Xj)^{xk,xi)'E{Aij^tigp)Ak,i,t{gp))\ ■ (5.32) 

i,j,k,UAN ) 


First note that, thanks to (|5.3ip . whenever i I and j ^ I, E (^Aij^t{9)Ak^i^t{s)) = 
^ ^ ^)) = 0, so that we can suppose that the indices {i,j,k,l) in 

()5.32p are such that i = I or j = 1. 
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We first treat the case i = 1: here, we can further assume that k = j, since if k j, one 
has that E{AijAk^i) = Ei(Ak^iE(Aij\fr,r + j) = 0. Hence, for fixed i,j, 


E ^i^i,jA9p)Aj,i,ti9p)) ^ E ^{^i,jA9pA)^'^{^jM9pA)A 

p ^ 1 P ^ 1 

p ^ 1 


E Y. '^{ 9 p{Tj,t,Ti,tf)\ , 


I P ^ 1 


\P ^ 1 




where S is defined in (j3.22j) . Thanks to Proposition 13.101 and Proposition I2.6I the last 
quantity is bounded by a constant that is uniform in i,j. Consequently, the contribution 
of the case i = I, j = k to (I5.32h can be bounded above by 


C 


(: 




\ \^n\ jjeAjv 


E ^{Xi,Xj) 


ki 


‘‘N 


|A.iv| IjeAjv 


E d{xi,Xj) 


-2a 


(5.33) 


When a < |, the latter quantity is bounded above by I ) = ^, by (I3.3jl . When 
a > i, this contribution is of order ( 

We now treat the case where j = I in (I5.32h . Proceeding as before and using again (I3.3jl , 
one easily sees that the contribution of this case to (I5.32p can be bounded above by 


C 




(5.34) 


which is uniformly bounded in N, in both cases a < | and a > |, by definition of a^v. One 
can conclude that (I5.32p is uniformly bounded in N. 

We now treat the case of the smallest contribution to E ^ j which corresponds 

to the case where all the four couples of indices are equal (see the case (G) in Figured]): 
h = h = h = h = (bi)- This case boils down to estimating the following quantity 


E ^iXi,XjA E ^{^i,jA9r?A,j,t{9sA) = E ^iXi,XjAE(^\\Aij^t\\tq) 

IAnI ijeAjv r,s^l |A7v| ijeAjv 

(5.35) 

E d{x.,x,)-^-. 

IAatI ijeAjv 


Using again we see that if 0 ^ a < |, the last quantity is bounded above by 

4 

Gj^j^ = 7 ^- If a = |, we obtain an upper-bound of order If | < a < i, the upper- 
bound is of order ^ ^ and in the case | < a < 1, one has an upper-bound of order 
^^ In any case, this term is bounded in N. From all this we can conclude 
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that 


sup E I 

t<T 




(5.36) 


'N,t 

We now turn to the second term in (15.271) . Setting [g,^t]{'r,x) = [g,^t]i{T,x) = 
J g{T,6,uj,x)^t{dS, do)) (recall (|2.2I) 1 and defining Aj := [xj,Xj+i[c S, j e A^r, one obtains 




~ Ta— i ^{xi,Xj) [g,it\ {Ti,t,Xj) - f ^{xi,x) [g,^t] {Ti^t,x)dx 

’ I^AfUeAjv \ ieAjv ) 

E L i^i^i^^j)[9,(t]{fi,t,Xj)-^{xi,x)[g,^t]{fi,t,x)) dx, 

I^A'I i,jeAj^9Aj 


CiN 

I^A'I ijeA 


E / [ 9 ,^t]{Ti^t,Xj){'i>{xi,Xj)-'^{xi,x)) dx 

lA I E / '^ixi^^)i[9,Ct]{fi^t,Xj)-[g,^t]{fi^t,x)) dx, 
IAatI jj.A.WA,' 


qn 


•“ '^nI(9) + '^n]:(9) 


'N,t 

r(3) 


Concerning j)^j, setting Uij{g) := [ 5 ,^*] {Ti^t,Xj) {^{xi,Xj) - ^'(xi,x)) dx, we have 

'^N^t _ ) “ 71 [4 E E '^iljli9r)ui2j2i9r)ui3j3igs)ui^j^{gs)\ , 

' ^ / LA. /V ?■ 1 7 . ^ \ r. .<? > 1 / 


“ 9 / |Ajv| \r,s ^ 1 

ilvJ4 


/ 






E E E 


1^2 

2 \ 


IAatI ij,A:,ZeAjv [ \p ^ 1 


I P ^ 1 


Using M, it is easy to see that |'I’(xi,Xj) - 'I'(xj,x)| dx ^ +i . Moreover, 

3 r”j| 

E ^{[9p,^t]{n,t,Xjf) ^ E f 9p{fi,t,0,oj,Xj)^Ct{d9,du:)\, 

p>i p >1 W / 


p 

E 


(/ 


5 - 3- 


^t{de, dta)] , 
< } 


which is uniformly bounded in i, j and N (recall Proposition 13.10]) . Consequently, for some 
constant C > 0 , 


(p ^ 1 

Thus, one obtains that 


E E ^m(5p)' E ^k,i{9pf\\ ^- 




sup E I 

t < T 




(p ^ 1 






E 


jY2(a-l) 






‘‘N 


iV4(l-a) ’ 


2 
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which is bounded, uniformly in N, by definition of oat. Turning to and setting 
■= = /a,. ^{xi,x){[g,^t] - [c/,6] {n,t,x)) (recall the defi¬ 

nition of the form U in Proposition 13.11]) . we have 


Hence, using the boundedness of the linear form U in (I3.28p 


E(ii 4 ‘iirj 


6 c- 




^ EfJ (1 + sup -x| / ^{xi^,x)dx\, 


IAtvI l=l 

iivj'4 


l xiAi 


^ c 




E Yli L '^ix^,x)dx\'E,Yl{l + \0^,tr^+ \uiX") ^ 

,...,14 l=l / l=l 




E ri E f >tte„i)dx<cA 

h,...,i4 1=1 jeAff 




so that sup^r ^ <j r E ^ j <+ 0 O. 

The same estimate about the last term /Cat,* can be proven along the same lines as we 
have done for the first term lN,t- We leave the proof to the reader. Proposition 13.121 is 
proven. □ 


5.5. Proof of Proposition 14.IL 

Definition 5.1. Define the linear forms on the space of functions {9,lj,x,6,uj,x) 
g{6,uj,x,6,u;,x) on (XxYxS)^; 


Vg{6,uj,9,u},x) '.= g{9,uj,x, 9, ui,x), 
Q9-= 9- Vg. 


(5.37) 

(5.38) 


If we fix (0, oj,9,uj, x) and we see g a.s a, function of the second spatial variable x only, we 
have to think of Vg as the first (constant) term in the Taylor decomposition of x fi((x) 
around x, Qg being the remainder. 

Proof of Proposition \4-l\ We decompose 7?Af,o 111.91) and PLn,d 112.7p into the sum of terms 
of different scalings (see Remark 15.21 belowf: for all test function (9,uj,x) /(^)Ci;,x) 


(f?7V,o ,/) = 7^ E {f{^i,0,<^UXi)-[f,^o]{Xi)) 

I^Af| i€AN 


“^(jA^ E [f^Co]{xi)-J [/,6](a;)dxj := (r/JJJ,,+ /), 


(5.39) 
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where [/, ^o] (x) = dcj) and for every test functions (r, f) i--> g(T,f) 

g) = E ^{Xi,Xj){g{Tifi,Tjfi) -[g,(o]{Ti^o,Xj)) 

I^Afl ijeAjv 

+ E (ta^ E ^ixi,Xj)[g,(o]{Ti,o,Xj)- fji>{xi,x)[g,Co]{Tifl,x)dx], 

l^ivUeAjvVl^ivI jeA^ -'S / 

■■={'^n1^9) + {'H^Z^9), (5.40) 

where (recall ([22])) [g,^o]{T,x) = [g,^o]i{9,<^,x,x) = 9{x,0,uj,x)^o{d9, dw). 

Remark 5.2. In (|5.39p and (|5.40l) . and 'H^q capture the fluctuations of the i.i.d. 
particles and disorder at t - 0 (and henee, should seale as for large N with Gaussian 

limits in the case on = yN). On the other hand, the processes and ^)vo capture 

the variations of the system w.r.t. the spatial variables (for a fixed randomness). Hence, 
(2) (2) 

the scaling of ryj^o o,nd 'Hy^-'o governed by the regularity of the weight At and the test 
functions. This scaling is in any case different from the Gaussian scaling and yields to 
deterministic limits as N ^ cx>. 

We first make the observation that in (I5.39p . ryyyg converges to 0 as oo, regardless 
of the value of a e [0,1). Indeed, since / in (j5.39p is at least of class (recall p3.9p ). one 
^ C II V/ Iloo which goes to 0 as ^ oo in both cases a < | and a > |. 
Secondly, using in (Ib.dOh the decomposition of Definition 15.11 we write 

[nf ,,, g) = , Vg) + , Qg). (5.41) 

The intuition for p5.4ip is that the only (possibly) nontrivial contribution to TLfiQ only 
comes from Vg, the remainder Qg being already sufficiently regular in {x, x) to compensate 
for the singularity of the kernel {x,x) df{x,x) in (I5.4nh : regardless of the value of 
ae [0,1), 



^ 0, as N ^oo. 

Indeed, since g is at least of class (recall p3.9ll i. one can write 

[Qg,f.o]i^,^^,x,x) = / (g{e,uj,x,e,0,x) - g{e,uj,x,e,u;,x))^oide, do), 

JXxY 

= f f dug{9,uj,x,9,0,u)duf,o{d6, do), 

^XxY Jx 


(5.42) 


d{x,x)i—^— f [dug,io\i{9,uj,x,u)du\, 
\a{x^x) Jx j 


\~Q(Q ^x ^x) 
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where, by assumption on g, the function Q is bounded and of class C^. Writing Qi{xi,Xj) = 



), {n 

a,' 

Qg'j can be written as 





(wS?b as) = 

ON 

I^Afl 

E 1 

(K E Kfe. 

V l^ivl jsAj^ 

Xj) [QffKo] 

{dgo, 


"j) 





- [Q9,^o] {9i,o, 

OJi^ Xi 


)■ 


= 

UN 

l^wl 

E 1 

ieAjv 

(lA 1 E d{xi,: 

V l^ivl jeA^ 

Tj)^-^QiiXi 

,Xj)- 


xx)^~"Qi{- 

Vi,x)dx^, 

= 

ttN 

l^wl 

E 1 

ieAjv 

[ E L {d{xi, 

VieAjv -^^9 

xA~°'Qi{x, 

i,Xj) 

-d{xi, 

x)^-"Qiixi, 

x)) dxj. 

= 

UN 

I^Af| 

E 1 

i£A]\i 

( E f. 

VjeAjv -^^9 

Xj) (^d{xi,x 

\ l-a 
?/ 

-d{xi, 

dxj 


+ 

“A 

I^aI 

■?i 


A~"iQiixi., 

,Xj)- 

■ Qii^Xi^ 

x)) dxj. 


= 


+ K'’o.9) 





(5.43) 


where Aj ;= [xj,Xj+i). We treat the two terms in (|5.43p separately. First of all, 


Kb 9) 


^ II d„. 


I^nI 


2 (E X 

jeAjv V jeAiv - 


\d{xi,Xj)^ "I dx 


} 


A long but easy calculation shows that the term within the parentheses is of order 
uniformly in i e A^r (it is the speed of convergence of the Riemann sum associated to 
towards its integral), so that ’ d) 

^ 00 . As far as the second term in (j5.43p is concerned, one has 


1 1 -a 


X \X^ 
to 0 as A" 


^ (7^, which, in any case in a, goes 


(<’ 0 . 9 ) 





that is of order ^ -^n^oo 0 too. This proves (|5.42l) . 

We specify now our analysis to cases a < | and a > |. We first treat the case a < |. 
First observe that in this case the whole process (^^ 0 ’ ^d)) converges to 

0 as A -5- 00 . Indeed, since one has trivially [Vg,^o] (r, x) = [ 5 ,^ 0 ] we have 


(Ww. P 9 ) 


1 

I^nI 


E 

ieAjv 


[ 5 , 6 ] {Tifl,Xi) 



d){Xi,Xj) 



which by Lemma 13.41 goes to 0 as A 00 . From this and the remarks made above, 
everything boils down to the identification of the limit of the process ^ i- 

To do so, we identify the limit of '■= + v , g'j for any scalars u and v 

and any test functions / and g. By a density argument, it suffices to identify the limit 
for test functions g of the type g{9,uj,x,9,uj,x) = gi{9,Lo)g2{9,uj)h{x,x). In such a case. 
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(recall that oat = \/N and IAtvI = 2A") 


{w5Jo,5) = i E 92{Oj,o,ujj){^. E gi{Oi,o,^^i)h{xi,Xj)'^{xi,Xj)], (5.45) 


where 


92(0,u}) := 92(0,u) - J §20,oj)^o{de, dw). 


(5.46) 


Intuitively, the term within brackets in (I5.45P should be, for large N, close to H(xj) where 
E{x) := - gi{9,uj)h{x,x)'d>{x,x)^o{dO, duj) dx. (5-47) 


Hence, one can decompose 77^0 into 

= i E 92{Oj,o,ujj)E{xj) 

VN j^An 

+ ^ E 52(%,o,Wj)(“(xj)-^ E 9i{di,o,^^i)h{xi,Xj)^{xi,Xj)] 

V N \ ieAjv / 

Consequently, Xw may be written as 

Xn = u , /) + u ,9) + v , 5 ), 

= y= E +u\/2g2(6'i,o,Wj)H(xi)j+ 5(| (5.49) 

Let us admit for a moment that 77^ g converges in to 0, as -> 00 . Then, Proposi¬ 
tion 14.11 for a < 4 is a consequence of Lyapounov Central Limit Theorem (see [8] , Theo¬ 
rem 27.3, page 362) applied to (|5.49l) and the convergence as A" -> 00 of 

E + v\/2g2{9ifi,uij)E{xi) 

j^An \Vv2 

to u‘^Cjj{f,f) + 2uvCri,n{f:9) + C'}i{g,g) (recall (j4.2l) and (l5.47p L We are now left with 
proving that 77^ g in ()5.48p goes to 0 in as A" ^ 00 . Indeed, using that {9ifl,iOi)i^AN 
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are i.i.d. and that £(52(^1,01 <^*)) = 0 for all i e Ajv, we obtain 

^ / iV 

“771^ E E(52(6'j,o,Wj)^)E(5i(6'fc,o,Wfc))^(a;fc,a:i)^'(xfc,a:i)H(xj) 

gfceAjv 

E E(52(6'i,o,Wj)^)E(5fi(6»j-o,Wi))E(cii(6»fc,o,Wfc))^(a;i,a:i)^'(xj,Xi)h(xfc,a:i)^'(xfc,Xi) 

ij-fceAjv 

= E(92(«i,„,o.,)")[^ Z 

ieAjv 

-E ( 51 ( 01 , 0 ,^ 1 )) E h{xk,Xi)'i/{xk,Xi)E{xi) 

i.fceAjv 

+ E( 5 i( 0 i,o,wi))^— ^ E /i(a:j,a:j)^(a7i,a7i)/i(a:fc,Xi)^'(xfc,Xi) 

As N ^ 00 ^ this quantity goes to 

7V^^((^^-0 ’J ^{x)‘^dx-AE{gi{6i^o,uJi))‘^ J J h{x,x)^{x,x)dx^ dx 

+ 2 E (5i(0i,o,a;i))^ fill h(x, x)'I'(x, x) dxj dx 
= 2 y" H(x)^dx-4 J H(x)^dx + 2 J H(x)^dx = 0. 

which is the desired result. This concludes the proof of Proposition 14. 1 1 in the case a < |. 

It remains to treat the case a > | (where aw = A7^““): as already proven in the case a < 

I, the quantity ZzeA^ (/(^*,o,Wi,Xi) - [/,^o] {xp) converges in law, so that = 

^ 1 ^ EifAjv {fi^i,o,‘xii,Xi) - [f,Co] (xi)) converges in law to 0 as iV ^ 00 . Since it is also 
the case for j^^o > /)> whole process 5^,0 converges in law to 0. It remains to identify 
the limit for T-Ln,o- 0 ’ d) ~ '^AO^S')’ '''^^ere J^]. was defined in (|5.27l) . In 

particular, we see from (I5.36|) and the definition of aw in (II.IUI) that vanishes to 

0 as Ai ^ 00 when a > ^. By the remarks made at the beginning of this proof, it only 
remains to identify the limit of Vg'^ in (I5.41h . An application of Lemma 13.41 to 

(j5.44p clearly shows that , Vg^ converges as A^ ^ 00 to the quantity dehned in (14.3p . 

This concludes the proof of Proposition 14.11 □ 

5.6. Proof of Proposition [4T6l The proof of this result is standard and uses techniques 
developed by Kunita in [2l] and used by Mitoma [32] and Jourdain and Meleard [23] . 

Since the following is mostly an adaptation of the approach of [23] to our case, we sketch 
the main lines of proof and refer to the mentioned references for technical details. 

Since 0 i-> c{9,uj) and P are supposed to be bounded as well as their derivatives up to 
order 3P+9, the function 6 i-> v{t,6,uj) is regular and bounded (with bounded derivatives), 
uniformly in t e [0,T]. According to [21], Theorem 4.4, p.277, the flow Xg^tix, f)o 
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defines a diffeomorphism. By [24j . for any differential operator D on (M™)^ of order 

smaller that 3P + 8 


sup sup (|Z)Xs,t(r,f)n < +0O, r > 0. 

(T,f)£XxYxS 0 ^ s sg t ^ T 


(5.50) 


Using this estimate and backward Ito’s formula ( |24] . Theorem 1.1, p.256), it is possible 
to prove (see [23], p.761) that for all g e C^, for all 0 ^ s ^ t ^ T, for all (r,f). 


Uit,s)g{T,f) - g{T,f) = £ ^y\U{t,r)g){T,f)dr. 


(5.51) 


The next step is to prove that (|5.51l) is also valid in the space C' 3 (p^ 2 ) • relies on the 
following lemma (we refer to |23] . Lemma 3.11 for a proof of this result): 

Lemma 5.3. Under the assumptions of Section \2.Sl the operator is continuous from 


^s'p+S ^3(P+2) 




a 


-0,2t 

3(P+2) 


^ C|| 5 f|Lo,. , [ 0 ,i] 




^o,2f. ^ ^II5^ II(7*^’^ 1 ^ ^1 ’ s^t € [0,r] . 

^3(P+2) 


For any j ^ 3P + 8 , the operator U{t, s) is a linear operator from Cj’^ to Cj’’' such that 

II U{t, s)g |Lo.i ^ C II <711^0,0 , 0 ^ s ^ ^ T, 

I U (t, s)g - U (t, s')g ||^o.i, ^ C || s' ||^o,o \/ s' - s, 


Note in particular that there exists C > 0 such that for any differential operator D of 
order smaller than 3(P + 2), || Z1 [TT, Vs\ - D [FT, nt] ||^ ^ C |t - s|, for any s,t e [0,T]. 
The change in the parameter l in the spaces Cj’'' in Lemma 15.31 comes from the fact 
that {0,uj) i-> c{9,u}) is possibly unbounded in to: there exists C > 0 such that for all 
differential operators Di and D 2 , with sum of orders smaller than 3(P + 2), = 

^ assumption on c. As a consequence of Lemma 15.31 


we 


know that if <7 e C'^p+g, s 1 -^ ^f‘^\U{t,s)g) is continuous in and hence that 


Ai), 


Jq ^j^\u{t, s)g) ds makes sense as a Riemann integral in C'^’^p^^y In particular, we 


0,2t 


obtain, for every g e 


3P-I-9 


U{t,s)g - g = £ ^y\U{t,r)g)dr, in Cg) 


2 t 

3(P+2)- 


(5.52) 


We are now ready to establish the representation (|4.18p : let S a solution to (j4.17|l in 


C([0,r],w^:7;^:7). Since 


2(P-l-2)' 


“JpIJ ^ have, for all g e Wf^-p^^y 


{£t,g) = £{Ss,^j^^g) ds + i: {Rs , g) ds. 


(5.53) 


Both relations (|5.52p and (|5.53p are in particular true for every g e Cgp^g, since C^’^p^, 2 ) 
^ 3 (P-i- 2 ) (I3.17P and since by assumption K>m and i> n + 2i, recall (I2.26P and (l2.27p L 
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Combining (I5.52|) and (I5.53|) , one obtains for g e C^p^g 


{£t,9) = + 

= £{i^y^r£s+Rs, 
= f^'{{^Pr£s+Rs, 


U{t,s)g-£ ^^^\U{t,r)g)dr'^ ds 

U{t, s)g) ds - JJ {i^s^^'^r£s + Rs , ^y\U{t,r)g)) dr ds, 
U{t,s)g) ds-£ fj {i^y'^r£s + Rs,^y\U{t,r)g)) dsdr, 
U{t,s)g^ ds- l^£r, ^yHu{t,r)g)'^ dr, (5.54) 


where we used again (I5.53|) in (15.54^ . Since C^p^g is dense in C' 2 (p+ 2 )’ 
fgU{t,s)*Rsds holds in C'° 2 '(p+ 2 )- Since C' 2 (p^ 2 ) dense in ^ 2 {P%y 
for (I4.17P in C([0,r], Proposition 14.61 is proven. 


the identity £t = 
uniqueness holds 
□ 
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